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1  INTRODUCTION 


BACKGROUND 

Nondestructive  testing  (NDT)  devices  are  being  widely  used  to 
evaluate  the  load-carrying  capability  of  pavements  for  air  carrier  and 
highway  pavements.  Evaluation  procedures  have  been  developed  using  vari¬ 
ous  types  of  NDT  devices  for  these  pavements.  Phase  1^  of  this  study 
evaluated  commercially  available  NDT  devices  for  use  on  light  aircraft 
pavements  (design  gross  loadings  less  than  30,000  lb).  Phase  II,  re¬ 
ported  herein,  is  to  develop  a  methodology  for  evaluation  of  light  air¬ 
craft  pavements  based  upon  multilayered  elastic  models  and  limiting 
stress/strain  criteria. 

PURPOSE 

The  purpose  of  Phase  II  of  this  study  is  to  develop  an  evaluation 
procedure  based  on  a  multilayered  elastic  procedure  for  evaluating  pave¬ 
ments  that  support  aircraft  with  gross  weights  of  less  than  30,000  lb. 

The  evaluation  will  determine  the  allowable  gross  aircraft  load  for  a 
given  number  of  operations. 

SCOPE 

This  study  will  utilize  data  only  from  nondestructive  testing 
devices  similar  to  the  Model  2008  Road  Rater.  Although  concepts  for 

the  model  are  general  and  would  probably  apply  to  any  of  the  other 
devices  evaluated  in  the  Phase  I  study,  this  device  alone  was  used  in 
development  of  the  evaluation  methodology. 


DEVELOPMENT  OF  THE  EVALUATION  METHODOLOGY 

REVIEW  OF  NDT  EVALUATION  PROCEDURES 

Since  a  number  of  NDT  pavement  evaluation  procedures  have  been 
developed,  a  cursory  review  will  be  made  to  outline  the  confines  for 

which  this  methodology  was  restricted. 

2 

Green  and  Hall  reported  a  procedure  that  uses  the  U.  S.  Array 
Engineer  Waterways  Experiment  Station  (WES)  l6-kip  vibrator.  This  pro¬ 
cedure  uses  a  dynamic  load  sweep  at  a  constant  frequency  of  15  Hz.  Only 
the  center  deflection  is  measured  and  the  resultant  parameter  (dynamic 
stiffness  modulus  (DSM))  is  computed  as  the  inverse  slope  of  the  upper 
third  of  the  load  deflection  relationship.  The  DSM  is  correlated  di¬ 
rectly  with  allowable  single-wheel  load.  The  procedure  utilizes  the 
U.  S.  Army  Corps  of  Engineer  (CE)  California  Bearing  Ratio  (CBR)  design 

for  flexible  pavements  and  the  Westergaard  theory  for  rigid  pavements. 

3 

Weiss  reported  a  layered  elastic  evaluation  procedure  using  re¬ 
sults  from  the  WES  16-kip  vibrator.  This  procedure  is  used  to  predict 
the  subgrade  modulus.  Since  the  Model  2008  Road  Rater  produces  a  i*000- 
lb  peak  load  (8000-lb  peak  to  peak)  as  compared  with  15,000-lb  peak  for 
the  WES  16-kip  vibrator,  this  procedure  could  not  be  readily  adapted. 

Yang  reported  a  procedure  that  uses  the  WES  l6-kip  vibrator  and 
a  frequency  sweep  where  the  dynamic  load  is  held  constant  and  the  fre¬ 
quency  is  varied  between  approximately  5  and  50  Hz.  This  procedure 
predicts  the  subgrade  modulus. 

The  frequency  sweep  test  requires  approximately  10  times  the 
amount  of  time  to  conduct  as  a  load  sweep.  WES  research  indicates  that 
the  data  derived  from  a  frequency  sweep  do  not  warrant  the  extra  time 
required  to  collect.  Pavement  material  properties  are  frequency- 
dependent.  Therefore,  another  variable  is  introduced  to  the  evaluation 
method  that  must  be  accounted  for. 

Treybig  et  al.^  reported  a  layered  elastic  evaluation /overlay  de¬ 
sign  procedure  using  the  Dynaflect  testing  device.  In  this  procedure, 
the  modulus  values  of  the  upper  layers  are  assumed.  Samples  are  taken 
of  the  granular  base  and  subbase  and  subgrade  materials,  which  are 
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tested  in  the  laboratory  using  a  dynamic  load  ti  i. axial  toot,  to  determine 
the  resilient  modulus.  By  use  of  the  deflection  from  the  number  one 
sensor  of  the  Dynaflect  and  assumed  or  laboratory  modulus  values,  the 
modulus  of  the  subgrade  is  determined  through  a  series  of  nomographs. 

This  modulus  is  compared  with  the  laboratory  results  of  resilient  modu¬ 
lus  versus  deviator  stress.  Through  a  series  of  approximations  of 
stress  and  modulus  that  parallels  the  laboratory  relationship,  the  design 
modulus  is  predicted.  The  basic  concepts  of  the  Treybig  procedure  are 
similar  to  the  approach  being  reported  except  that  a  goal  in  developing 
this  procedure  is  not  to  require  destructive  sampling  and  laboratory 
testing. 

Koole^  reported  a  layered  elastic  evaluation/overlay  design  pro¬ 
cedure  using  the  Falling  Weight  Deflectometer .  The  deflection  at  the 
center  of  the  applied  load  and  the  ratio  of  that  deflection  to  a  deflec¬ 
tion  some  distance  away  from  the  applied  load  were  used.  The  pavement 
was  characterized  by  a  three-layer  system.  The  subgrade  modulus  and 
either  the  asphaltic  concrete  (AC)  modulus  or  the  AC  thickness  are  deter¬ 
mined  in  this  procedure.  The  modulus  values  for  the  AC  layer  (if  the 
thickness  is  to  be  determined)  and  layer  2  are  determined  from  laboratory 
results  or  from  construction  records.  The  effective  thickness  of  the 
surface  layer  is  determined.  From  these  values,  the  layered  elastic 
theory  is  used  to  predict  allowable  or  overlay  requirements. 

•7 

Anani  reported  results  using  the  Model  1*00  Road  Rater.  The 
Model  boo  Road  Rater,  described  in  detail  in  Phase  1  of  this  study, 
applies  a  peak-to-peak  dynamic  load  of  approximately  720  lb  to  the 
pavement  through  two  4-  by  7-in.  pads  spaced  6  in.  apart.  Four  sensors 
are  used  to  monitor  deflections.  One  is  located  between  the  pads  while 
the  other  three  are  spaced  at  1-ft  intervals. 

Anani  used  the  BISAR  computer  program  developed  by  the  Shell  Oil 
Company  with  successive  approximation  techniques  to  determine  the  elastic 
moduli  of  the  pavement  layers.  Using  up  to  a  four-layer  system,  the 
procedure  approximated  the  E  values  by  the  following  equation: 

E(I)  =  E(I)  x  £MlLi_AilLRRD(I) 
new  old  2 


where 


E(I) 

new 

E(I)oia 

RRD(I) 


A(I) 


modulus  value  for  layer  I 

previous  assumed  modulus  value  for  layer  I 

measured  deflection  corresponding  to  layer  I,  i.e., 
deflection  No.  It  is  dependent  on  layer  ^  and  deflec¬ 
tion  No.  3  is  dependent  on  layer  3,  etc. 

deflection  from  BISAR  associated  with  E(I)  . 
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Other  procedures  have  been  developed  by  Sharpe  and  Ullidtz 
that  predict  modulus  values  for  pavement  systems  using  the  layered 
theory  and  nondestructive  testing  devices.  These  procedures  use  two 
deflection  measurements  from  the  deflection  basin. 

As  shown  above,  a  number  of  researchers  used  nondestructive 
testing  and  the  layered  elastic  theory  to  evaluate  and  design  overlays 
for  pavements.  The  procedure  developed  in  this  study  will  incorporate 
some  of  the  stronger  points  of  those  procedures  as  well  as  add  some  new 
approaches. 

EXPERIMENTAL  DATA 


The  data  used  in  developing  this  procedure  were  collected  on  the 
Pennsylvania  Transportation  Research  Facility  (PTRF)  and  on  selected 
pavements  at  the  WES.  The  PTRF  is  a  one-mile  track  located  at  Pennsylvania 
State  University.  Fifteen  items  were  tested  with  the  Model  2008  Road 
Rater  in  June  1978.  From  1  to  20  tests  were  conducted  on  each  item. 

The  deflections  used  in  this  study  will  be  the  average  of  those  tests. 

These  averages  are  given  in  Table  1  for  the  7000-  and  5000-lb  force 
levels.  The  pavement  structures  of  the  15  items  are  given  in  Table  2. 

Also  shown  in  Table  2  are  the  number  of  equivalent  l8-kip  axle  loads 
(EAL)  that  were  applied  to  the  pavement  items  prior  to  testing.  The 
surface  conditions  of  the  items  were  good  except  for  items  3,  A,  and  C. 
These  items  exhibited  some  surface  distress  such  as  longitudinal  and 
alligator  cracking. 

Two  test  pits  were  excavated  in  the  PTRF  to  obtain  samples  of  the 
granular  subbase  material  and  the  subgrade.  Undisturbed  block  samples 
and  bag  (disturbed)  samples  were  taken  of  the  subgrade.  Bag  samples 
were  also  taken  of  the  subbase.  Resilient  modulus  tests  were  conducted 
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on  these  samples  in  accordance  with  the  procedures  outlined  in  Appendix 
A.  Ave rage  results  of  the  subgrade  tests  are  presented  in  Figure  1. 

It  should  be  noted  that  undisturbed  samples  gave  higher  modulus  values 
than  disturbed  samples.  The  granular  subbase  results  are  presented  in 
Figure  2. 

RELATIONSHIP  OF  LAYERED  ELASTIC 
THEORY  TO  MEASURED  DEFLECTIONS 

GENERAL 

The  first  assumption  in  developing  this  procedure  is  that  dynamic 
deflections  correspond  to  those  deflections  predicted  from  the  layered 
elastic  theory.  To  validate  this  assumption  experimental  data  were  com¬ 
pared  with  results  from  two  computer  programs.  The  first  is  the  Shell 
BISAR  computer  program,  based  on  the  layered  elastic  theory,  which  re¬ 
lates  stress  and  strain  in  each  layer  to  a  load  applied  at  the  surface 
of  a  pavement.  The  other  program,  entitled  CHEVIT,10  gives  a  nonlinear 
approach  to  the  solution  of  the  modulus  of  the  lower  layers  by  using 
laboratory  stress-modulus  relationships  for  granular  and  subgrade  layers. 
Those  layers  are  divided  into  sublayers  for  which  the  stress  is  calcu¬ 
lated  initially,  and  from  the  laboratory  stress-modulus  relationship  a 
new  modulus  is  computed.  With  this  modulus,  the  program  again  computes 
the  stress.  The  program  iterates  until  a  solution  is  obtained  for 
modulus  and  stress. 

LINEAR  ANALYSIS 

To  determine  the  applicability  of  the  deflection  basin  to  the 
layered  elastic  analysis,  the  BISAR  computer  program  was  initially  used. 
The  modulus  values  for  the  AC  surface  layers  were  obtained  from  the  16- 
Hz  relationship11  presented  in  Figure  3.  An  evaluation  of  data  collected 
with  the  WES  l6-kip  vibrator  on  specially  constructed  temperature  sections 
at  the  WES  was  made  using  the  l6-Hz  modulus  versus  temperature  relation¬ 
ships.  The  test  sections  were  small  (approximately  20  by  20  ft)  and 
consisted  of  varying  thickness  of  AC  over  a  lean  clay  (CL)  subgrade. 

The  design  of  these  sections  considered  the  thickness  of  AC  to  be  the 
only  variable.  Tests  were  conducted  over  a  wide  range  of  temperatures. 
Modulus  values  for  the  AC  as  determined  from  Figure  3  were  used  in  the 
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Figure  3.  Modulus-temperature  relationships  for 
AC  surfacing  (l  psi  =  703  kg/m2;  1°F  =  -17°C) 

(after  Kingham  and  Kallas11) 

BISAR  program  to  predict  deflections.  Good  correlation  was  obtained  for 
the  temperature  test  section  between  the  deflection  ratios  (deflection 
at  a  given  temperature  to  deflection  at  70°F)  obtained  from  the  experi¬ 
mental  results  and  the  results  from  the  BISAR.  Therefore,  the  relation¬ 
ships  presented  in  Figure  3  were  selected  for  determination  of  the  AC 
modulus  values. 

Modulus  values  for  the  pavement  layers  other  than  the  AC  surface 
layers  of  the  PTRF  were  estimated  from  construction  data  and  based  on 
laboratory  results  that  were  collected  during  the  initial  construction 
(not  the  laboratory  results  presented  in  this  report).  Poor  agreement 
was  obtained  between  the  BISAR  deflections  and  the  measured  deflections 
for  the  PTRF  sections.  Therefore,  it  was  concluded  that  it  is  extremely 
difficult  to  estimate  correct  modulus  values  for  pavement  layers  that 
have  been  subjected  to  extensive  traffic  and  environmental  effects. 

NONLINEAR  ANALYSIS 

An  analysis  was  then  made  using  the  CHEVIT  nonlinear  program. 
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The  granular  subbase  layer  was  characterized  by  the  laboratory  relation¬ 
ships  shown  in  Figure  2.  The  subgrade  resilient  modulus  relationships 
in  Figure  U  represent  the  averages  of  the  undisturbed  sample  tests. 
Figure  5  shows  a  typical  pavement  section  (PTRF  item  1A).  The  initial 
estimate  for  the  modulus  of  the  nonlinear  layers,  as  well  as  the  modulus 
for  the  linear  layer,  is  given.  Pavement  sections  for  all  of  the  PTRF 
items  are  shown  in  Table  2. 

The  summation  of  the  strains  in  the  bottom  layer  to  infinity  by 
the  layered  elastic  model  tends  to  give  larger  deflections  than  the 
measured  values.  To  compensate  for  this  effect,  a  rigid  layer  was 
placed  in  this  system  model  at  a  depth  of  20  ft  below  the  surface. 

Figure  6  shows  a  comparison  of  predicted  deflections  with  and  without 
the  rigid  boundary  to  the  measured  deflections.  Note  that  the  basins 
predicted  from  CHEV1T  using  the  rigid  layer  agree  better  with  the  Model 
2008  Road  Rater  basins  than  those  predicted  without  the  rigid  layer. 


Figure  U.  Relationship  of  subgrade  resilient 
modulus  and  deviator  stress  used  in  CHEVIT 
(l  psi  =  703  kg/m2) 


(LINEAR) 


Figure  5.  PTRF  item  1A  pavement  section  breakdown 
for  CHEVIT  input  (l  in.  =  2.5*4  cm;  1  psi  =  703  kg/m  ) 
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Figure  6.  Comparison  of  layered-elastic  deflection  basins 
with  and  without  rigid  layer  to  measured  basins  (l  in.  = 
2.5^  cm:  1  mil  =25.1  microns;  1  lbf  =  U.418  N) 


! 

Basins  are  shovm  for  each  of  the  15  test  items  of  the  PTRF  in 
Figures  7  through  11.  Tests  with  both  5000-  and  7000-lb  loads  were  con¬ 
ducted  and  compared  with  CHEVIT  results.  Note  the  good  agreement  in  the 
basin  data  for  those  sections  without  the  lime  and  cement-treated  base 
courses.  The  modulus  values  for  the  lime  and  cement-treated  layers  were 
taken  from  construction  records.  The  heavy  traffic  on  the  PTRF  at  the 
time  of  test  apparently  resulted  in  cracking  of  these  stabilized  layers. 
Therefore,  the  stabilized  layers  had  much  lower  moduli  values  at  the 
time  of  the  test  than  those  values  used  in  the  CHEVIT  analysis.  This 
fact  explains  the  variance  in  measured  and  predicted  deflections  in 
items  3,  A,  C,  E,  F,  and  G  shown  in  Figures  8,  10,  and  11.  From  the 
results,  it  appears  that  the  deflection  basin  is  a  measurement  that  can 
be  modeled  with  a  layered  elastic  theory  and,  therefore,  used  to  predict 
the  strength  parameters  of  the  pavement  layers. 

EFFECTS  OF  STATIC  PRELOAD 

The  effects  of  the  static  load  applied  to  the  pavement  surface  as 
a  preload  with  the  Model  2008  Road  Rater  were  analyzed  using  the  CHEVIT 
program.  Figure  12  illustrates  loading  versus  time  for  the  Model  2008 
Road  Rater  for  the  5000-  and  7000-lb  tests.  The  sinusoidal  loading  ap¬ 
plies  a  minimum  and  a  maximum  force  to  the  pavement  surface.  The  magni¬ 
tudes  of  these  forces  range  from  500  to  6500  lb  for  the  5000-lb  peak-to- 
peak  force  level  and  from  500  to  7500  lb  for  the  7000-lb  peak-to-peak 
force.  The  CHEVIT  program  was  run  with  forces  of  6500  and  1500  lb,  and 
the  differences  (diff)  in  predicted  deflections  from  each  run  were 
calculated  to  model  the  sinusoidal  loading  of  the  Road  Rater.  The  re¬ 
sults  were  compared  to  a  run  where  the  force  was  5000  lb.  Table  3  shows 
the  results  of  these  calculations  for  PTRF  items  1A,  IB,  and  2.  A  maximum 
difference  of  6.3  percent  occurred  at  the  5000-lb  force  level  at  12  in. 
from  the  center  of  the  load;  however,  the  percent  difference  was  prac¬ 
tically  negligible  for  most  comparisons.  This  analysis  indicates  that 
the  effect  of  the  static  load  for  computer  modeling  of  the  Road  Rater 
results  is  negligible  particularly  when  the  vibrator  is  operated  near 
the  maximum  output.  Therefore,  the  static  load  will  be  neglected  in  the 
determination  of  the  layer  modulus  values. 


DISTANCE  FROM  CENTER.  IN  DISTANCE  FROM  CENTER.  IN  DISTANCE  FROM  CENTER.  IN 


SUM  *NO 1 1D31 J30 


(l  in.  =  2.5U  cm;  1  mil  =  25. ^ 


LOAD.  KIPS 


FTRF  1A 


Diff  7000  6.382  5.1*86  4.042  1.7C 

( Continued ) 


4.440  N:  1  in.  -  cm 


Table  3  (Concluded) 


Load 

Deflections,  mils 
Cited  Distances,  : 

,  at 
in. 

lb 

0 

9 

12 

18 

21* 

36 

1*8 

60 

6.507 

5.59*1 

1+.91+0 

3.778 

2.901* 

1.790 

1.107 

0.772 

Percent 

Di f ference 

2.0 

2.0 

2.0  2.0 

PTRF  IB 

1.7 

1.0 

0.0 

-0.5 

7.128 

5.829 

5.011* 

3.598 

2.61+1 

1.589 

0.970 

0.691+ 

^■im 

1.81*9 

1.511 

1.299 

0.922 

0.660 

0.372 

0.217 

0.153 

Diff  5000 

5.279 

1* .  318 

3.715 

2.676 

1.981 

1.217 

0.753 

0.51+1 

5000 

5  -  6lL 

1* .  588 

3.9*+l+ 

2.820 

2.057 

1.221 

0.739 

0.528 

Percent 

Difference 

6.3 

6.3 

6.2 

5-1+ 

3.8 

0.3 

-1.9 

-2.1* 

7500 

8.120 

6.61*1* 

5.718 

1+.113 

3.030 

1.837 

1.127 

0.806 

500 

0.6U6 

0.1*55 

0.323 

0.229 

0.126 

0.072 

0.050 

Diff  7000 

7.1*7U 

6.115 

5.263 

3.790 

2.801 

1.711 

1.055 

0.756 

7000 

7.625 

6.237 

5.366 

3.856 

2.835 

1.713 

1.01*8 

0.750 

Percent 

Difference 

2.02  - 

2.0 

2.0  1.7 

PTRF  2 

1.2 

0.1 

0.6 

-0.8 

6500 

6.071 

5.116 

1+.1+92 

3.1+05 

2.613 

1.618 

1.001 

0.696 

1500 

1.1+71 

1.238 

1.086 

0.817 

0.620 

0.371+ 

0.226 

0.157 

Diff  5000 

1* .  600 

3.878 

3.1+06 

2.588 

1.993 

1.21*1+ 

0.775 

0.539 

5000 

U.706 

3.963 

3.1+77 

2.630 

2.013 

1.21*0 

O.76U 

0.532 

Percent 

Difference 

2.3 

2.2 

2.1 

1.6 

1.0 

-0.3 

-1.1* 

-1.3 

7500 

6.983 

5.887 

5.171 

3.925 

3.018 

1.871+ 

1.161 

0.807 

0.502 

0.1+23 

0.279 

0.211 

0.126 

0.052 

Diff  7000 

6.1*81 

5.1+61* 

1* .  800 

3.61*6 

2.807 

1.71+8 

1.086 

0.755 

7000 

6.527 

5-501 

1* .  831 

3.665 

2.815 

1.71+6 

1.080 

0.751 

Percent 

Difference 

0.7 

0.7 

0.6 

0. 52 

0.3 

-0.1 

-0.6 

-0.5 

DETERMINATION  OF  LAYER  MODULUS  VALUES  FROM  BASIN  DATA 
DEVELOPMENT 

The  deflection  basin  produced  by  applying  a  load  to  the  pavement 


with  the  Model  2008  Road  Rater  gives  four  input  parameters  to  the  system 
analysis  that  can  be  used  to  derive  the  strength  parameters  of  the  pave¬ 
ment  layers.  A  program  called  CHKVDEF  was  developed  to  determine  a  set 
of  modulus  values  that  provide  the  best  fit  between  a  measured  deflection 
basin  and  a  computed  deflection  basin  when  given  an  initial  estimate  of 
the  modulus  values,  a  range  of  modulus  values,  and  a  set  of  measured  de¬ 
flections  . 

Consider  the  pavement  system  where: 

a.  The  modulus  is  unknown  for  a  number  of  layers  (NL). 

b.  The  deflection  due  to  plate  load  is  measured  at  a  number  of 
deflection  (ND)  locations. 

c.  ND  is  greater  than  NL  . 

The  objective  is  to  determine  the  set  of  E's  that  will  minimize 
the  error  between  the  computed  deflection  A  and  the  measured  deflection 
RRD  .  To  accomplish  the  objective,  a  relationship  was  developed  for  the 
deflection  at  a  point  j  as  a  function  of  the  unknown  E's  ,  i.e., 

Aj  =  f(E1,  Eg. • .Enl) 

then  the  error  at  a  position  where  the  deflection  was  measured  is 


RRD.  -  A  =  RRD,  -  f(E?,  E2,  .  .  .  E^ ) 


This  expression  is  then  squared  and  summed  with  respect  to  each  measured 
deflection 


To  minimize  the  error  with  respect  to  an  unknown  E  ,  the  partial  deriva¬ 
tive  of  the  error  function  is  taken  with  respect  to  the  E  .  By  taking  a 
derivative  with  respect  to  each  unknown  E  ,  then  a  set  of  NL  equations 
is  obtained  that  can  be  solved  giving  the  set  of  E's  for  the  minimum 


error  between  the  measured  basin  and  the  computed  basin. 

First,  a  set  of  E  values  is  assumed  and  the  deflection  A°  is 

tJ 

computed  corresponding  to  the  measured  deflection  RRD  .  Each  unknown 

J 
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E  is  varied  individually  and  a  now  set  of  deflections  is  computed  for 
each  variation.  lining  the  two  computed  deflections  and  the  two  values 
of  each  E  ,  a  function  is  determined  for  each  deflection.  For  example, 
let 


E*  =  log10  E 


Then  the  deflection  at  location  1  is  given  as  a  function  of  E  ,  i.e.. 


A1  =  A11  +  S11E^1 


where 


A°  A1 

A1  ‘  A1 

Sl1  =  E£°  -  EH* 

A11  =  A1°  -  S11E? 

Ei.°  =  logjQ  of  first  assumed  value  of  E^ 

E£^  =  log10  ° f  ®i  after  the  variation 

A°  =  computed  deflection  at  position  1  for  E° 

A^  =  computed  deflection  at  position  1  for  E^ 

Likewise,  functions  are  determined  for  each  deflection  and  each  unknown 
E  ,  resulting  in  j  =  1  to  ND  and  i  =  1  to  NL  .  Then 


A.  =  A . .  +  S..EJI. 
J  Ji  Ji  i 


To  write  an  expression  for  A.  as  a  function  of  all  E's  ,  the  fol¬ 


J 


lowing  is  used 


Aj  =  A°  +  (changes  in  Aj  due  to  changes  in  the  E's) 

Consider  when  the  modulus  of  layer  changes  from  E°  to  ,  the  change 

-i  _  J.  J. 

in  Aj  would  be  Sj^(E£^  -  ES.^ )  . 


Thus 


NL 

A,  =  +  £  S,.(Ei.  -  El°) 

j  j  ji  i  a 
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The  value  of  A^  can  be  expressed  in  terms  of  any  of  the  unknown  E's  , 
i.e.,  ,  as 

AJ  =  A.1NL  +  ‘"’.INLAND 

The  expression  for  A .  now  becomes 

J 


IM-Lj 

\j  =  AJNL  +  SjNLEiNL  +  ^1Sji(E*'i  “ 


The  error  squared  for  the  jth  position  is  (RRD  -  A  )  or 

0  o 


ERROR*:  = 
J 


(  r  o  NL  )2 

b  -  [aj»l  *  Vlm.  *  f=1sJi(Eli  -  EtP] 


The  summation  of  the  error  for  all  readings  is 
ND  ND  (  r 

T  ERROR.  =  y  <RRD,  -  A,.,t  +  S  E£ 

7=1  J  j=i{  J  L  jNL  JNL 


A +  S  E£ 
jNL  jNL  . 


ML 

°  +  y  S..(E*.  -  EH°m 

NL  i-lJ1  1  lJf 


If  a  weight  term  for  each  reading  is  to  be  applied,  then  the  expres¬ 

sion  becomes 

ND  .  ND  /  ,  NL  A 

^[Wj  "  (ERR0R)]  =  i ^Wj |RRDj  "  [AjNL  +  SjNLEilNL  +  ^=1Sji(EAi  “  Efci)]|j 

Taking  the  partial  with  respect  to  each  E  and  setting  the  partial  equal 
to  zero,  the  following  is  obtained: 

ND  f  r  NTj  i  s 

° '  £>wj  b  ■  [V  +  V&  +  L  sJi(EL  -  EI?>] 

If  the  equations  derived  are  put  in  the  form 

[B]  {E}  =  {0} 

the  {C}  terms  are  the  constant  part  of  the  equation.  For  k  =  1  to  NL 


ND 

y  s„w. 

RRD,  - 

('a.,. 

+  S  „  ES,° 

NL 

_  y 

f=i Jk  J 

J 

\  JNL 

JNL  NL 

i-: 

for  k  =  1 

to  NL 

and 

i  =  1  to  NL 

is 

B,  .  =  j  8..W.S,. 
ki  4*=1  jk  J  ji 


2k 


r 


If  the  weight  term  is  chosen  to  he 
developing  the  equation  from 

ERROR j  = 


Wj  =  RRDj 


RRD,  -  A. 

_J _ 1 

RRD 

J 


the  result  is  the  same  as 


which  is  a  percent  type  error.  The  solution  of  the  equation  is  the  set 
of  E's  that  minimizes  the  percent  error.  The  efficiency  of  the  proce¬ 
dure  will  depend  on  how  well  the  functions  represent  the  actual  relation¬ 
ship  between  the  computed  deflection  and  the  E's  . 

It  appears  that  as  long  as  the  final  E  values  are  within  the 
initial  input  limits,  the  A  =  fClog^^  E)  is  a  good  representation  of 
the  relationship. 

A  computer  program  named  CHEVDEF,  consisting  of  the  procedure 
described  above,  was  used  in  developing  the  pavement  evaluation  proce¬ 
dure  reported  herein.  CHEVDEF  uses  the  CHEVRON  layered  elastic  program 
as  a  subroutine  to  compute  surface  deflections.  A  flowchart,  input 
format,  example  input,  example  output,  and  a  listing  of  CHEVDEF  are 
presented  in  Appendix  B. 

The  limitations  of  this  approach  are  that  the  layered  elastic 
theory  assumes  a  uniform  pressure  applied  to  the  surface  of  the  pavement. 
With  the  Model  2008  Road  Rater,  the  load  is  applied  through  a  rigid  cir¬ 
cular  plate  with  the  center  deflection  measured  on  top  of  that  plate. 
Therefore,  a  difference  does  exist  in  the  measured  center  deflection  and 

a  deflection  computed  from  layer  elastic  procedures  at  the  center  of  the 

12 

load  area,  as  illustrated  by  Fossberg  in  Figure  13.  Note  that  the 
elastic  layer  solution  and  field  data  coincide  at  approximately  three- 
fourths  of  the  radius  of  the  plate.  In  order  to  determine  the  optimum 
spacing  for  the  deflection  measurements,  computations  such  as  shown  in 
Figures  14  and  15  were  made  with  the  CHEVDEF  program  for  flexible  and 
rigid  pavements,  respectively.  Spacing  distances  of  0,  4-5,  6,  and  9 
inches  were  used  for  the  computations.  Varying  these  distances  caused 
little  change  in  the  subgrade  modulus  for  either  the  flexible  or  the 
rigid  pavement.  There  was  variation  in  the  surface  modulus  on  the 
flexible  pavement.  Based  on  the  temperature-frequency  relationship 
presented  in  Figure  3,  the  modulus  of  the  asphalt  layer  should  be 
520,000  psi.  Thus,  from  Figure  15,  a  distance  between  4.5  and  6  in. 
appears  to  produce  the  best  results. 
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RADIAL  OFFSET.  IN. 


Figure  13-  Patterns  of  resilient 
deflections  under  a  n . -diameter 

plate  (plate  pressure  =  h  psi  )  di¬ 
rectly  on  the  subgrade  (l  in.  = 

2.5^1  etn)  (after  Fossberg^) 


SENSITIVITY  ANALYSIS 

To  evaluate  tbe  accuracy  of  the  CHEVDEF ,  an  analysis  was  made  on 
a  thin  and  a  thick  pavement,  section.  The  CHEVRON  program  was  used  to 
calculate  deflections  for  the  two  pavement  sections.  These  deflections 
were  used  as  measured  deflections  for  this  analysis  so  that  an  error 
associated  with  field  measurements  would  be  eliminated.  Tables  h  and  5 
present  the  results.  The  modulus  values  listed  as  correct  values  in 
these  tables  were  input  to  the  CHEVRON  program  to  give  the  deflections 
that  were  used  as  measured  deflections  in  CHEVDEF.  In  the  first  case, 
the  initial  estimate  for  modulus  was  higher  than  the  correct  values. 

The  initial  estimate  was  lower  than  the  correct  values  for  the  second 


Reproducibility  of  the  CHEVDEF  Program  for  a  Thin  Section 
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703  kg/m  ;  1  in.  =  2.5^  cm;  1  mil  =  25-U  microns. 


case.  In  tin;  third  cane,  t.ho  initial  cr.t  imalc  was  a  1  lernat  ed  between 
low  and  high  values.  The  program  reasonably  reproduced  all  modulus 
values  for  the  thick  section  but  varied  from  the  correct  values  for  the 
thin  upper  layers  for  the  thin  section.  The  values  for  the  the  two 
bottom  layers  were  very  close  even  for  the  thin  section. 

To  determine  the  reproducibility  for  different  spacings  of  deflec¬ 
tions,  the  data  for  the  thick  section  were  selected.  Table  6  lists  the 

results  of  this  analysis.  No  significant  difference  occurred  due  to  the 

location  of  the  deflection  with  respect  to  the  final  value  of  modulus 
for  up  to  four  layers. 

Another  check  was  made  on  a  typical  section  with  a  stiff  (stabi¬ 
lized)  layer.  Table  7  summarizes  the  results.  Again  the  program  pro¬ 
duced  a  good  approximation  of  the  correct  values  for  the  four-layer  sys¬ 
tem. 

To  determine  the  sensitivity  of  the  CHEVDEF  program  to  the 

Poisson's  ratio  of  each  layer,  an  analysis  was  made  using  the  thick  sec¬ 

tion.  The  Poisson's  ratio  was  varied  between  0.2  and  0.5  for  each 
layer.  Figure  16  shows  the  change  in  modulus  values  for  these  varia¬ 
tions  in  the  Poisson's  ratio.  There  is  little  effect  on  the  predicted 
modulus  values  when  the  Poisson's  ratios  for  layer  2  and  layer  3  are 
varied.  The  variations  are  large  in  layer  1  and  layer  2  when  the 
Poisson's  ratios  of  layer  1  and  layer  are  varied.  In  all  cases, 
there  is  little  variation  in  the  modulus  of  the  lower  layers. 

A  Poisson's  ratio  of  0.35  will  be  assigned  to  all  layers  above 
the  subgrade,  and  a  ratio  of  O.ii  will  be  assigned  to  the  subgrade. 

PREDICTION  OF  LAYER  MODULUS  VALUES 

The  CHEVDEF  program  was  used  to  predict  modulus  values  for  the  15 
FTRF  test  items.  Table  8  presents  the  results.  Layer  1  includes  all 
the  AC  material,  and  the  modulus  value  for  this  layer  war.  taken  from 
the  l6-Hz  relationship  in  Figure  3.  The  pavement  temperature  was  com¬ 
puted  from  the  surface  temperature  plus  the  previous  five-day  mean . ^  ~ 

Basins  predicted  from  the  program  were  not  plotted  against  the 
measured  basins,  but  the  small  differences  can  be  seen  in  the  column 
entitled  "Absolute  Sum  of  Differences  in  Deflections."  In  each  case, 
this  is  the  sum  of  the  differences  in  four  deflections. 
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Table  7 

Reproducibility  of  the  CHEVDEF  Prog 
for  a  Section  with  a  Stiff  Laver 


Figure  16.  Effect  of  varying  Poisson's  ratios 
vl’  v2’  v3’  and  vl4  )  on  predicted  modulus  values 
(l  ksi  =  6.89  MPa) 
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Table  8 

Predicted  Modulus  Values  for  PTRF  Items  free  CHEVPEF 
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The  stress  produced  in  the  subgrade  by  the  Model  2008  Road  Rater 
may  not  be  the  stress  associated  with  the  design  aircraft  load.  To  ac¬ 
count.  for  thin,  two  loads  were  analyzed  (r>000  and  7000  lb).  The  modulus 
values  and  stresses  (vortical,  longitudinal,  and  rudial)  at  the  top  of 
the  subgrade  were  calculated  for  each  load  level.  With  these  values,  a 
relationship  between  the  modulus  and  the  deviator  stress  can  be  developed 
for  stress-dependent  subgrades.  Table  8  also  shows  the  slope  of  the 
line  and  the  intercept,  as  well  as  the  bulk  stress  at  the  top  of  the 
subbase  layer.  There  appears  to  be  no  relationship  to  the  laboratory 
results  for  these  values.  The  values  of  E  for  the  granular  layer  are 
values  that  satisfy  the  model  deflection  basin.  The  granular  material 
without  a  binder  has  very  few  load  transfer  properties  through  bending. 
Therefore,  the  model  predicted  lower  modulus  values.  Associated  with 
the  low  modulus  value  is  a  high  vertical  strain  in  that  layer,  which  is 
considered  a  problem  in  using  the  layered  elastic  theory  to  model 
granular  materials. 

DETERMINATION  OF  SUBGRADE  MODULUS  FOR  EVALUATION 

Results  from  the  CHEVDEF  program  give  the  relationship  for  the 
deviator  stress  and  the  modulus  for  the  subgrade  materials  in  the  form 


E  =  SoD  +  I 


where 

E  =  subgrade  modulus 
S  =  slope 

=  deviator  stress 
I  =  intercept 

It  is  known  from  Hooke's  Law  and  by  definition  for  resilient  modulus  that 

E  =  — 

e 

where  e  represents  the  strain.  By  substituting  and  rearranging 


E  =  SEe  +  I 


1  -  Se 


If  a  limiting  vertical  compressive  strain  is  selected,  the  modulus  of 
the  subgrade  can  be  calculated.  For  comparisons  used  in  this  study, 
a  value  of  0.0006  in. /in.  was  selected.  This  value  was  taken  from  an 
average  line  drawn  on  Figure  IT  and  represents  500,000  repetitions  or 
25,000  arrival /departures  per  year  for  20  years.  The  equation  for  the 
average  line  is  given  as 

c  =  ANB 
v 

where 

ev  =  limiting  vertical  compressive  strain  in  subgrade,  in. /in. 

A  =  0.0063548 

N  =  number  of  strain  repetitions 
B  =  -0.17985 

The  modulus  values  in  Table  9  associated  with  that  strain  level  will  be 
used  for  the  evaluation  of  these  pavements. 

COMPARISON  OF  SUBGRADE  MODULUS  VALUES 

The  design  modulus  can  also  be  computed  from  laboratory 
resilient  modulus  test  results  as  outlined  in  Appendix  A.  The  best- 
fit  line  for  deviator  stresses  between  5  and  12  psi  was  calculated  and 
then  presented  in  Table  9-  A  comparison  of  the  mean  subgrade  modulus 
values  from  the  CHEVDEF  program  analysis  ( 1*4,201  psi)  to  the  laboratory 
results  (12,864  psi)  seems  very  reasonable. 

Another  comparison  was  made  to  the  subgrade  strength  parameters 

2 

derived  from  the  DSM  method  reported  by  Green  and  Hall.  The  PTRF  fa¬ 
cility  was  tested  with  the  WES  l6-kip  vibrator  at  the  same  time  that 
testing  was  conducted  with  the  Model  2008  Road  Rater.  Table  10  sum¬ 
marizes  results  of  the  DSM  evaluation.  The  subgrade  strength  factor  is 
the  equivalent  of  the  design  CBR  for  the  subgrade.  The  factor  is  calcu¬ 
lated  as  the  CBR  required  to  support  a  single-wheel  load  with  a  254-sq- 
in.  contact  area  for  1200  annual  departures.  The  load  on  the  wheel  is 
determined  by  a  correlation  with  the  measured  DSM.  The  relationship  of 
1500  times  CBR  was  used  to  compare  the  design  modulus  values  for  the  DSM 
method  and  for  the  layered  elastic  procedure  (Figure  l8).  Except  for 
two  items  (PTRF  items  4  and  B)  that  had  very  low  deflections  and  high 
DSM' s  due  to  stabilized  layers,  subgrade  strength  design  parameters 
approximate  each  other. 
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Deviator  Stress  -  Resilient  Modulus 
Relationships  from  Laboratory  Tests 


Sample 

Type 

Slope 

Intercept 

Design  E* 

No. 

Sample 

S 

I 

psi 

1-1 

Remolded 

-1279 

18,1*28 

10,1*21* 

1-2 

Remolded 

-1052 

15,632 

9,585 

II-l 

Undisturbed 

-ll8l 

20,811 

12,182 

11-2 

Undisturbed 

-  932 

22,581 

ll*  ,1*81* 

III-l 

Remolded 

-33**2 

1*5,998 

15,305 

III-2 

Remolded 

-10U9 

15,089 

9,26l 

III-3 

Undisturbed 

-1773 

38,808 

18,80)+ 

Mean 

=  12,861* 

Standard  Deviation 

=  3,512 

Note:  1 

p 

psi  =  703  kg/m  . 

*  Based 

on  strain  level  of 

0.0006 

in. /in. 

Table  10 

DSM  Evaluation  of  the  PTftF 


Temperature 

Allowable  Gross 

Corrected 

Subgrade 

Aircraft 

Load,  kips 

PTRF 

DSM 

Strength 

25l+  sq-in. 

127  sq-in. 

Item* 

kips/ in. 

Factor 

Single  Wheel** 

Single  Wheelt 

1 A 

1210 

6.95 

11.3 

86.8 

IB 

1008 

9-29 

92.7 

79-9 

2 

1331 

18.97 

122.1* 

103.8 

3 

813 

11-79 

71*. 7 

58.3 

1* 

2286 

26.1* 

210.3 

181.9 

5 

863 

9.58 

79.>* 

61.6 

6 

151*2 

15-9 

11*1.9 

116.8 

7 

125!+ 

15.2 

115.1* 

96.9 

8 

657 

9.05 

60.1* 

1*5-5 

A 

605 

11.1 

55-7 

1*0.2 

B 

163I* 

25-6 

150.3 

131.9 

C 

66l 

8.85 

60.8 

1*5-5 

E 

708 

6.1*9 

65.1 

1*8.2 

F 

859 

9.5 

79.0 

61.3 

G 

689 

9.15 

63.3 

1*7 . 5 

Note: 

1  kip  =  1*. 1*1*8 

kN;  1  kip/in. 

=  1.75  kN/cm;  1 

sq  in.  =  6.1*5 

sq  cm. 


*  See  Table  2  for  a  description  of  these  pavements. 

**  1,200  annual  departures,  20-year  life, 
t  25,000  annual  departures,  20-year  life. 


Figure  l8.  Comparison  of  subgrade  modulus 
by  the  DSM  and  layer  elastic  method  (l  ksi  = 
6.89  MPa) 

EVALUATION  OF  ALLOWABLE  AIRCRAFT  LOADS 


For  the  evaluation  of  the  allowable  aircraft  loads,  the  PAVEVAL 

3 

computer  program  reported  by  Weiss'  was  used.  For  flexible  pavements, 

this  program  calculates  the  vertical  compressive  strain  at  the  top  of 

the  subgrade  and  the  tensile  strain  at  the  bottom  of  the  AC  layer. 

These  strains  are  compared  to  the  limiting  strain  criteria  reported  by 

1*4 

Barker  and  Brabston.  For  rigid  pavements,  the  limiting  tensile  stress 
in  terms  of  the  number  of  load  (stress)  repetitions  from  Parker  et  al.1'’ 
is  considered  as 

_ R 

°RL  A  +  B  log  (cov) 

where 


o  =  limiting  value  of  tensile  stress,  psi 
nL 

R  =  flexural  strength,  psi 

A  =  0.58901 
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B  =  0.35^36 

cov  =  number  of  coverages  (The  number  of  coverages  is  determined 
by  dividing  the  number  of  aircraft  departures  by  the 
departure-to-coverage  ratio.  The  ratio  used  for  single¬ 
wheel  gears  was  7.9^  and  for  dual-wheel  gears  was  5.2.) 

Example  inputs  and  outputs  f'or  evaluation  of  light  aircraft  allowable 
loads  on  rigid  and  flexible  pavements  are  given  in  Appendix  C.  A  com¬ 
parison  of  the  allowable  gross  loads  for  the  PTRF  items  was  made  be¬ 
tween  the  DSM  method  and  the  layered  elastic  procedure  (Figure  19). 

For  comparison  purposes,  the  design  loads  were  calculated  for  a  single¬ 
wheel  aircraft  with  127-sq-in.  contact  area  and  25,000  annual  departures. 


ALLOWABLE  GROSS  LOAD,  KIPS 


0* - 1 - 7 - 1 - 1 - , 

0  25  50  75  100  12 

LAYERED  ELASTIC  THEORY 

Figure  19 .  Comparison  of  allowable  gross 
aircraft  loads  from  the  DSM  and  layered 
elastic  methods  (l  kip  =  b . I4U8  kN) 
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For  the  layered  elastic  evaluation,  tile  modulus  ot‘  the  AC  was  se¬ 
lected  as  770,000  psi,  which  corresponds  to  flu-  70°F  temperature  in 
Figure  3-  The  modulus  of  the  subbase  and  base  course  layers  (E2  in 
three-layer  items,  and  E2  and  E3  in  four-layer  items)  were  taken  from 
the  7000-lb  load  in  Table  8. 

The  allowable  loads  calculated  from  the  DSM  method  are  generally 
higher  than  those  calculated  using  PAVEVAL.  The  controlling  strain  is 
shown  in  Table  11  for  each  of  the  PTRF  items.  The  lower  allowable  loads 
for  the  layered  elastic  evaluations  are  probably  due  to  the  tensile 
strain  controlling  and  being  based  on  the  initiation  of  a  crack  in  the 
bottom  of  the  AC  rather  than  the  crack  propagating  through  the  AC  to  the 
surface  as  was  the  case  in  the  DSM  evaluation  (CBR  Design  System). 


Table  11 

Allowable  Gross  Aircraft  Loads 
from  PAVEVAL 


Allowable  Gross 

Controlling 

Strain 

PTRF 

Aircraft  Load 

Vertical 

Tensile 

Item 

kips 

Subgrade 

AC 

1A 

5U.7 

X 

IB 

bb.2 

X 

2 

52.6 

X 

3 

37.9 

X 

k 

92.6 

X 

5 

82.1 

X 

6 

77-9 

X 

7 

79-8 

X 

8 

hb.2 

X 

A 

6l.l 

X 

B 

120.0 

X 

C 

23.2 

X 

E 

27. h 

X 

F 

31.6 

X 

G 

31.6 

X 

Note:  1  kip  =  kN. 

A  limited  amount  of  testing  was  conducted  on  rigid  pavements  to 
verify  the  applicability  of  this  procedure.  Tests  with  the  Model  2008 
Road  Rater  were  conducted  on  a  Portland  cement  concrete  (PCC)  road  with 
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the  pavement  structure  as  shown  in  Figure  20.  Laboratory  resilient  modu¬ 
lus  tests  were  not  conducted  on  the  subbase  but  were  conducted  on  the 
subgrade  material  (Figure  21). 

Table  12  presents  the  CHEVDEF  program  results  for  this  rigid  pave¬ 
ment.  The  subgrude  modulus  selected  for  evaluation  for  ri,iT,id  pavements 
is  the  modulus  that  Ln  associated  with  a  5-psi  deviator  stress  and  is 
calculated  directly  from  the  slope  and  intercept  derived  from  the  devia¬ 
tor  stress-modulus  relationship  for  the  5000-  and  7000-lb  loads. 

Also  shown  in  Table  12  is  the  modulus  for  5-psi  deviator  stress 
from  the  laboratory  resilient  modulus  test.  CBR  tests  were  taken  to 
depths  of  30  in.  The  sample  taken  from  the  30-  to  U5-in.  depth  indicated 
a  higher  strength  material  (Figure  20).  Therefore,  the  predicted  sub¬ 
grade  modulus  seems  to  be  a  reasonable  estimate  of  the  subgrade  strength. 


NONDESTRUCTIVE  EVALUATION  PROCEDURE 


GENERAL 

The  procedure  outlined  in  this  section  is  for  both  rigid  and  flex¬ 
ible  pavements  and  is  based  on  a  layered  elastic  model  that  characterized 
multilayered  pavement  systems.  The  layer  strength  parameters  are  com¬ 
puted  from  field  in  situ  measurements.  The  strength  parameters  will  be 
input  into  an  evaluation  program  that  is  designed  to  handle  multiwheel 
aircraft  at  varying  traffic  levels.  The  output  will  be  computed  as  the 
allowable  load  for  a  20-year  design  life  pavement.  The  evaluation  will 
be  valid  for  conditions  existing  at  the  time  of  test  and  will  not  account 
for  changes  due  to  such  factors  as  frost  action  or  moisture  in  the  sub¬ 
grade.  These  factors  should  be  accounted  for  through  conventional  proce¬ 
dures  . 

NONDESTRUCTIVE  TESTING  EQUIPMENT 

The  need  for  a  device  with  variable  loading  characteristics  is 
required  to  describe  the  nonlinear  characteristics  of  the  subgrade  mate¬ 
rial.  Laboratory  resilient  modulus  tests  would  be  required  to  describe 
these  nonlinear  characteristics  when  the  description  cannot  be  derived 
from  the  NDT  data. 

The  NDT  device  must  output  a  minimum  of  three  deflections  of  which 
one  is  measured  near  the  applied  load  and  the  others  are  spaced  to  a  dis¬ 
tance  of  at  least  36  in.  from  the  applied  load.  The  number  of  deflec¬ 
tion  measurements  limits  the  number  of  variable  modulus  layers  that  can 
be  analyzed  with  this  procedure.  It  is  recommended  that  the  magnitude 
of  the  first  deflection  (that  deflection  measured  nearest  to  the  applied 
load)  be  at  least  0.2  mil. 

DATA  COLLECTION 

PAVEMENT  INFORMATION 

Before  evaluating  a  light  aircraft  pavement,  information  as  to 
pavement  types,  layer  thicknesses,  and  layer  types  must  be  derived  from 
construction  records  or  from  destructive  tests  (cores  or  test  pits). 

This  information  is  required  for  the  evaluation  of  allowable  load  as  well 
as  for  the  determination  of  NDT  locations.  If  test  pits  or  cores  are 
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r'.qui  rod ,  these  tests  :;h<>uJd  Ik-  per  I'urmed  after  the  NDT. 
a.vi  low  deflections  should  hi-  included  with  areas  selected  through  a 
conventional  sampling  procedure  to  identify  the  causes  of  these  unique 
IJDT  responses.  Rigid  pavement  cores  or  beams  should  be  tested  for 
flexural  strength. 

The  pavement  condition  should  be  surveyed  to  determine  the  areas 
and  types  of  distress.  It  is  not  within  the  scope  of  this  procedure  to 
outline  a  detailed  condition  analysis.  Since  material  properties  are 
affected  by  water,  frost,  and  primarily  environmental  conditions,  it  may 
be  necessary  to  reduce  the  strength  parameter  for  anticipated  freeze/thaw 
conditions,  alligator  cracking  in  flexible  pavements,  and  the  presence 
of  joints  and  cracks  in  rigid  pavements. 

TEMPERATURE  DATA 

Temperature  data  are  required  for  AC  pavements  to  evaluate  the 
pavement  properties  at  the  time  of  testing  and  the  ability  of  the  pave¬ 
ment  system  to  support  future  operations.  The  pavement  temperature 
during  the  tine  of  testing  is  determined  at  the  mid-depth  of  the  AC 

layer,  the  pavement  surface,  ambient  temperature,  and  the  previous  five- 

13 

day  mean  air  temperature  from  Figure  22.  The  pavement  surface  tempera¬ 
ture  should  be  measured  at  one-hour  intervals  during  the  period  of 
testing.  The  mean  air  temperature  can  be  obtained  from  the  nearest 
office  of  the  National  Oceanographic  and  Atmospheric  Administration. 

To  evaluate  the  pavement  system  for  future  operations,  the  average 
daily  maximum  temperature  and  average  daily  mean  temperature  is  needed 
for  the  hottest  months  and  the  spring  thaw  months.  These  data  can  also 
be  obtained  from  the  nearest  office  of  the  National  Oceanographic  and 
Atmospheric  Administration. 

NDT  TEST  DATA 

NDT  should  be  conducted  at  100- ft  intervals  alternating  to  either 
side  of  the  center  line  in  the  wheelpaths  on  flexible  pavements.  Rigid 
pavement  tests  should  be  conducted  in  the  slab  center  also  alternating 
to  either  side  of  the  feature  center  line.  A  minimum  of  five  tests 
should  be  conducted  in  each  pavement  type.  Parking  aprons  should  be 
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PAVEMENT  SURFACE  TEMPERATURE  PLUS  5-DAY  MEAN  AIR  TEMPERATURE, °F 


tested  on  a  200-ft  grid  system.  Testing  should  not  be  conducted  when 
the  pavement  or  subgrade  is  frozen. 

The  NOT  device  should  be  operated  at  the  frequency  that  produces 
t.he  bent  signal  response.  For  the  Model  POOR  Road  Rater,  this  frequency 
is  15  Hz.  Deflection  basin  data  should  be  collected  at  the  maximum  force 
output  of  the  device  and  at  a  force  level  of  50  to  75  percent  of  the  maxi¬ 
mum  force  output. 

DETERMINATION  OF  LAYER  MODULUS  VALUES 

The  CHEVDEF  program  provides  a  tool  with  which  the  modulus  values 
of  up  to  four  layers  can  be  predicted.  The  CHEVDEF  program  input  guide, 
typical  input,  and  program  listing  are  furnished  in  Appendix  B.  Sensi¬ 
tivity  studies  presented  earlier  in  this  report  showed  that  the  moduli 
of  the  subgrade  and  other  lower  layers  in  the  pavement  system  reproduce 
well  even  when  the  surface  moduli  may  differ.  Also  since  the  model 
assumes  a  uniform  pressure  and  the  Model  2008  Road  Rater  applies  a  rigid 
plate  to  the  pavement,  it  is  recommended  that  the  modulus  values  of  the 
surface  layer  be  assigned  as  described  below. 

For  flexible  pavements,  the  modulus  of  the  AC  should  be  determined 
from  Figure  3-  The  pavement  temperature  can  be  computed  from  the  surface 
temperature  plus  the  previous  five-day  mean  air  temperature  (Figure  22). 
The  depth  should  be  the  midpoint  of  the  AC  thickness. 

Table  13  gives  approximate  ranges  for  the  modulus  of  the  pavement 
materials.  These  ranges  represent  very  broad  limits  that  should  be 
compressed  for  the  CHEVDEF  program.  For  example,  if  a  pavement  with  a 
cement-stabilized  base  is  evaluated  and  the  traffic  history  indicates 
very  little  heavy  traffic  has  been  appLied,  a  range  for  the  stabilized 
base  should  be  selected  between  600,000  and  1,500,000  psi. 

NOT  results  should  be  analyzed  by  plotting  the  number  one  sensor 
deflection  versus  the  distance  along  the  feature.  Pavement  areas  should 
be  divided  according  to  pavement  types  and  pavement  deflections.  For 
each  area,  each  deflection  of  the  basin  should  be  averaged  to  determine 
the  mean  deflection  basin. 


2  — 

—  =  D0  ;  etc . 


where 

=  deflection  from  number  one  sensor 
n  =  number  of  tests  in  particular  items 
D  =  deflection  to  be  used  for  evaluation 

The  mean  deflection  values  for  each  pavement  area  should  be  input  to 
CHEVDEF  in  mils. 


Table  13 

Typical  Ranges  for  Modulus  Estimates 
and  Poisson's  Ratio  Values  for  Pavement  Layers 


Material 

Range  of  Modulus 
10^  psi 

Assigned  Value 
of  Poisson's  Ratio 

AC 

100  -  1000 

0.35 

PCC 

4000  -  6000 

0.15 

Untreated  base 

2  -  160 

0.35 

Treated  base 

8  -  2000 

0.35 

Subgrade 

2-37 

0.4 

2 

Note:  1  psi  =  703  kg/m  . 

*  Temperature-  and  frequency-dependent. 


In  analyzing  the  results  from  the  CHEVDEF  program,  it  is  impor¬ 
tant  to  check  the  predicted  modulus  for  a  layer  against  the  limit.  If 
the  modulus  does  hit  a  limit,  the  program  should  be  rerun  modifying  the 
limits  to  include  the  predicted  E  disregarding  boundary  conditions 
(see  Appendix  E,  Example  Output). 

The  highest  load  should  be  evaluated  first  with  the  CHEVDEF.  The 
values  for  the  modulus  for  the  upper  layers  obtained  from  this  run  may 
be  used  for  constant  values  when  running  the  lower  load  to  determine  the 
relationship  of  stress  versus  modulus  for  the  subgrade.  For  example,  if 
the  pavement  contains  three  layers,  consisting  of  AC,  granular  base,  and 
subgrade,  the  modulus  for  the  AC  and  granular  layer  as  determined  from 
the  initial  run  of  the  high  load  should  be  held  constant  at  the  values 
for  the  run  at  the  lower  load. 


Once  the  modulus  values  and  deviator  stress  are  obtained,  the 
following  equations  are  used  to  describe  the  nonlinearity  of  the  subgrade 


and 


S  = 


I 


where 

S  =  slope 

=  predicted  modulus  for  high  load  (7000  lb  for  2008  Road  Rater) 
E^  =  predicted  modulus  for  low  load  (5000  lb  for  2008  Road  Rater) 

=  deviator  stress  (vertical  -  radial  stress)  for  7000-lb  load 
=  deviator  stress  (vertical  -  radial  stress)  for  5000-lb  load 
I  =  intercept 

After  determining  the  S  and  I  for  the  equation  above,  calcu¬ 
late  the  limiting  strain  for  AC  pavement  from  Figure  17  based  on  the 
design  life  of  the  pavement,  and  compute  the  design  subgrade  modulus 
from  the  equation  given  below: 


E  = 


where 

E  =  design  subgrade  modulus,  psi 
c  =  limiting  strain 

For  rigid  pavements,  the  design  subgrade  modulus  will  be  selected 
at  a  deviator  stress  of  5  psi,  or 

E  =  5S  +  I 


DETERMINATION  OF  ALLOWABLE  AIRCRAFT  LOADS 

3 

The  PAVEVAL  program  presented  by  Weiss  will  be  used  to  predict 
allowable  aircraft  loads.  Example  inputs  for  both  rigid  and  flexible 
pavements  are  shown  in  Appendix  C.  Aircraft  characteristics1^  for  light 
aircraft  pavements  are: 


Type 

Gross 

Max /Load 
kips 

Departure  to 
Coverage  Ratio 

Contact 
Area 
sq  in. 

Wheel 

Spacing,  in 

Single 

20 

7.91* 

127 

— 

Dual 

30 

5-2 

75 

18 

Note:  1 

kip  =  U.UU8 

kN ;  1  sq  in.  = 

6.U5  sq  cm 

;  1  in.  = 

The  modulus  for  the  AC  surface  layer  should  be  determined  based  on  a  design 
pavement  temperature  for  input  to  Figure  3.  The  method  of  selecting  the 
design  pavement  temperature  for  this  evaluation  is  taken  from  Brabston 

n  #T  ■«  Q 

et  al.  Witczak  presents  a  relationship  between  pavement  temperature 
and  air  temperature  (Figure  23)  that  can  be  used  to  determine  the  design 
pavement  temperature  if  the  corresponding  design  air  temperature  is 
known.  For  this  design  procedure,  the  design  air  temperature  for  a  par¬ 
ticular  locale  is  determined  by  averaging  the  average  daily  maximum 
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temperature  and  the  average  daily  mean  temperature  for  the  design  month. 

Generally,  the  set  of  average  temperatures  will  be  necessary  only  for 

the  hottest  month  indicated  in  the  reporting  period. 

Thid  method  of  calculating  a  design  modulus  for  the  AC  layer  is 

conservative.  Other  procedures  have  been  reported  by  Barker  and  Brabston 

and  Koole.^  Modulus  values  of  those  layers  of  base  and  subbase  should 

be  selected  from  the  CHEVDEF  output  for  the  high-load  (7000- lb)  data. 

These  values  approximate  the  stress  of  the  light  aircraft  and  should  be 

representative  for  the  behavior  when  an  elastic  model  is  used. 

The  PAVEVAL  program  has  the  capability  of  calculating  AC  and  PCC 

3 

overlays.  These  procedures  are  available  in  the  Weiss  report  but  were 
not  evaluated  for  use  in  light  aircraft  pavement  design  as  a  part  of 
this  study.  The  computer  program  listed  in  Appendix  B  may  be  repro¬ 
duced  for  use  by  any  interested  party. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


This  study  has  resulted  in  the  development  of  an  evaluation  pro¬ 
cedure  for  light  aircraft  pavements  based  on  a  layered  elastic  model. 
Nondestructive  pavement  test  results  are  used  to  predict  the  layer 
strength  parameters  that  can  be  input  into  a  layered  elastic  model  to 
predict  the  allowable  load-carrying  capacities  of  both  rigid  and  flex¬ 
ible  pavements  containing  either  stabilized  or  nonstabilized  layers. 

The  deflection  basins  measured  from  the  NDT  device  at  two  force 
levels  are  used  for  input  to  a  system  that  predicts  the  nonlinear 
stress-dependent  behavior  of  the  subgrade  material.  Results  compare 
favorably  with  laboratory  resilient  modulus  tests. 

It  is  recommended  that  this  procedure  be  adopted  for  use  in  eval¬ 
uating  light  aircraft  pavements.  Further  study  should  be  conducted  to 
make  this  approach  applicable  to  air  carrier  airport  pavements.  The  use 
of  a  finite  element  code  for  the  modeling  of  a  rigid  plate  and  the  non¬ 
linear  behavior  of  the  granular  materials  should  be  included  in  future 
studies. 
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APPENDIX  A:  LABORATORY  PROCEDURE  FOR  DETERMINING 
THE  RESILIENT  MODULUS  OE  .'5IIHGRADK  SOUS* 


The  objective  of  this  test  procedure  is  to  determine  a  modulus 
value  for  subgrade  soils  by  means  of  resilient  triaxial  techniques.  The 
test  is  similar  to  a  standard  triaxial  compression  test,  the  primary  ex¬ 
ception  being  that  the  deviator  stress  is  applied  repetitively  and  at 
several  stress  levels.  This  procedure  allows  testing  of  soil  specimens 
in  a  repetitive  stress  state  similar  to  that  encountered  by  a  soil  in  a 
pavement  under  a  moving  wheel  load. 

DEFINITIONS 


The  following  symbols  and  terms  are  used  in  the  description  of 


this  procedure: 

a.  =  total  axial  stress. 

b.  =  total  radial  stress;  i.e.,  confining  pressure  in  the 

triaxial  test  chamber. 


— ’  °d  =  al  ~  °3  =  devia'tor  stress;  i.e.,  the  repeated  axial 
stress  in  this  procedure. 

d_.  e1  =  total  axial  strain  due  to  . 

£.  Mp  =  a^/e^  =  resilient  modulus. 

f.  0  =  +  2a  =  +  3°3  =  sum  °f  the  principal  stresses  in  the 

triaxial  state  of  stress. 


£.  °2//°3  =  principal  stress  ratio. 

h_.  Load  duration  =  time  interval  over  which  the  specimen  is 

subjected  to  a  deviator  stress. 

i..  Cycle  duration  =  time  interval  between  successive  applications 

of  a  deviator  stress. 


SPECIMENS 


Various  diameter  soil  specimens  may  be  used  in  this  test,  but  the 
recommended  specimen  diameter  is  2.5  to  3.0  in.  or  approximately  four 
times  maximum  aggregate  size.#  The  specimen  height  should  be  at  least 
twice  the  diameter.  Undisturbed  or  laboratory  molded  specimens  can  be 
used.  Procedures  for  obtaining  undisturbed  soil  specimens  are  given  in 

Hi 


*  This  procedure  is  taken  from  the  report  by  Barker  and  Brabston. 


Methods  for  laboratory 


Engineer  Manual  1110-2-1907*  "Uoil  Sampling. 

preparation  of  molded  specimens  and  for  back-pressure  saturation  of 

specimens,  if  required,  are  presented  in  EM  1110-2-1906,  "Laboratory 
20 

Soils  Testing." 

EQUIPMENT 

TRIAXIAL  TEST  CELL 

A  triaxial  cell  suitable  for  use  in  resilience  testing  of  soils 
is  sb'  *n  in  Figure  A-l.  This  equipment  is  similar  to  most  standard  cells 
with  che  exceptions  of  being  somewhat  larger  to  facilitate  the  internally 
mounted  load  and  deformation  measuring  equipment  and  having  additional 
outlets  for  the  electrical  leads  from  the  measuring  devices.  For  the 
type  of  equipment  shown,  air  or  nitrogen  is  used  as  the  cell  fluid. 

The  external  loading  source  may  be  any  device  capable  of  pro¬ 
viding  a  variable  load  of  fixed  cycle  and  load  duration,  ranging  from 
simple  cam-and-svitch  control  of  static  weights  or  air  pistons  to  a 
closed-loop  electrohydraulic  system.  A  load  duration  of  0.2  sec  and  a 
cycle  duration  of  3  sec  have  been  found  to  be  satisfactory  for  most  ap¬ 
plications.  A  square-wave  load  form  is  recommended. 

DEFORMATION  MEASURING  EQUIPMENT 

The  deformation  measuring  equipment  consists  of  linear  variable 
differential  transducers  (LVDT's)  attached  to  the  soil  specimen  by  a 
pair  of  clamps.  Two  LVDT's  are  used  for  the  measurement  of  axial  deforma 
tion.  The  clamps  and  LVDT's  are  shown  in  position  on  a  soil  specimen  in 
Figure  A-l.  Details  of  the  clamps  are  shown  in  Figure  A-2.  Load  is  meas 
ured  by  placing  a  load  cell  between  the  specimen  cap  and  the  loading 
piston  as  shown  in  Figure  A-l. 

Use  of  the  type  of  measuring  equipment  described  above  offers 
several  advantages: 

a.  It  is  not  necessary  to  reference  deformations  to  the  equipment 
which  deforms  during  loading. 

b.  The  effect  of  end-cap  restraint  on  soil  response  is  virtually 
eliminated. 
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£.  Any  effects  ui'  piston  friction  arc  eliminated  by  measuring 
loads  inside  the  triaxial  cell. 

In  addition  to  the  measuring  devices,  it  is  also  necessary  to 
maintain  suitable  recording  equipment.  It  is  desirable  to  have  simul¬ 
taneous  recording  of  load  and  deformation.  The  number  of  recording 
channels  can  be  reduced  by  wiring  the  leads  from  the  LVDT's  so  that  only 
the  average  signal  from  each  pair  is  recorded.  The  introduction  of 
switching  and  balancing  emits  permits  use  of  a  single-chamber  recorder. 
However,  this  will  not  permit  simultaneous  recording. 

ADDITIONAL  EQUIPMENT 

In  addition  to  the  equipment  described  above,  the  following  items 
are  also  used: 

a.  A  10-  to  30-ton-capacity  loading  machine. 

b.  Calipers,  a  micrometer  gage,  and  a  steel  rule  (calibrated  to 
0.01  in. ) . 

c.  Rubber  membranes,  0.01  to  0.025  in.  thick. 

d.  Rubber  o-rings. 

£.  A  vacuum  source  with  a  bubble  chamber  and  regulator. 

f,  A  back-pressure  chamber  with  pressure  transducers. 

£.  A  membrane  stretcher. 

h.  Porous  stones. 

PREPARATION  OF  SPECIMENS  AND  PLACEMENT  IN  TRIAXIAL  CELL 

The  following  procedure  should  be  followed  in  preparing  and 
placing  specimens: 

20 

a.  In  accordance  with  procedures  specified  in  EM  1110-2-1906, 
prepare  the  specimen  and  place  it  on  the  base-plate  complete 
with  porous  stones,  cap,  and  base  and  equipped  with  a  rubber 
membrane  secured  with  0-rings.  Check  for  leakage.  If  back¬ 
pressure  saturation  is  anticipated  for  cohesive  soils,  pro¬ 
cedures  indicated  in  Appendix  X  to  EM  1110-2-1906  for  the 
Q-type  triaxial  tests  should  be  followed.  For  purely  non- 
cohesive  soils,  it  will  be  necessary  to  maintain  the  vacuum 
during  placement  of  the  LVDT's.  The  specimen  is  now  ready 
to  receive  the  LVDT's. 

b.  Extend  the  lower  LVDT  clamp  and  slide  it  carefully  down  over 
the  specimen  to  approximately  the  lower  quarter  point  of  the 
specimen. 

£.  Repeat  this  step  for  the  upper  clamp,  placing  it  at  the  upper 
quarter  point.  Insure  that  both  clamps  lie  in  horizontal 
planes . 


d.  Connect  the  LVDT's  to  the  recording  unit,  and  balance  the  re¬ 
cording  bridges.  This  step  will  require  recorder  ad,' usrnents 
and  adjustment  of  the  LVDT  steins.  When  a  recording  bridge 
baiancf  has  been  obtained,  determine  (to  the  nearest  0.01  in.) 
the  vertical  spacing  between  the  I.VPT  clamps  and  record  this 
value . 

e_.  Place  the  triaxiai  chamber  in  position.  Set  the  load  cell  in 
place  on  the  specimen. 

£.  Place  the  cover  plate  on  the  chamber.  Insert  the  loading 
piston,  and  obtain  a  firm  connection  with  the  load  cell. 

g.  Tighten  the  tie  rods  firmly. 

h.  Slide  the  assembled  apparatus  into  position  under  the  axial 
loading  device.  Bring  the  loading  device  to  a  position  in 
which  it  nearly  contacts  the  loading  piston. 

i_.  If  the  specimen  is  to  be  back-pressure  saturated,  proceed  in 
accordance  with  EM  1110-2-1906. 

j_.  After  saturation  has  been  completed,  rebalance  the  recorder 
bridge  to  the  load  cell  and  LVDT's. 

RESILIENCE  TESTING  OF  COHESIVE  SOILS 

The  resilient  properties  of  cohesive  soils  are  only  slightly  af¬ 
fected  by  the  magnitude  of  the  confining  pressure  .  For  most  appli¬ 
cations,  this  effect  can  be  disregarded.  When  back-pressure  saturation 
is  not  used,  the  confining  pressure  used  should  approximate  the  expected 
in  situ  horizontal  stresses,  which  will  generally  be  on  the  order  of  1 
to  5  psi.  A  chamber  pressure  of  2  psi  is  a  reasonable  value  for  most 
testing.  If  back-pressure  saturation  is  used,  the  chamber  pressure  will 
depend  on  the  required  saturation  pressure. 

Resilient  properties  are  highly  dependent  on  the  magnitude  of  the 
deviator  stress  .  It  is  therefore  necessary  to  conduct  the  tests  for 
a  range  in  deviator  stress  values.  The  following  procedure  should  be 
followed: 

a. .  If  back-pressure  saturation  is  not  used,  connect  the  chamber 

pressure  supply  line  and  apply  the  confining  pressure  (equal 
to  the  chamber  pressure).  If  back-pressure  saturation  is 
used,  the  chamber  pressure  will  already  have  been  established. 

b.  Rebalance  the  recording  bridges  for  the  LVDT's,  and  balance 
the  load  cell  recording  bridge. 

£.  Begin  the  test  by  applying  1000  repetitions  of  a  deviator 
stress  of  not  more  than  one-half  the  unconfined  compressive 
strength. 
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d_.  Decrease  the  deviator  load  to  the  lowest  value  to  he  used. 

Apply  200  repetitions  of  load,  recording  the  recovered  verti¬ 
cal  deformation  at  or  near  the  last  repetition. 

e.  Increase  the  deviator  load,  recording  deformations  as  in 
Step  d.  Repeat  over  the  range  of  deviator  stresses  to  be 
used.  It  is  recommended  that  3,  5,  7,  12,  15,  17,  20,  22, 
and  25  j*-_i  be  used. 

f.  At  the  completion  of  the  loading,  reduce  the  chamber  pressure 
to  zero.  Remove  the  chamber  LVDT's  and  load  cell.  Use  the 
entire  specimen  for  the  purpose  of  determining  the  moisture 
content . 

The  results  of  the  resilience  tests  can  be  presented  graphically 
as  shown  in  Figure  A-3  for  the  resilient  modulus  and  in  the  form  of  a 
summary  table  such  as  Table  A-l. 


DEVIATOR  STRESS  ad  ,  PSI 


Figure  A-3.  Presentation  of  results  of  resilience 
tests  on  cohesive  soils  (l  psi  =  703  kg/in  ) 
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RESILIENCE  TESTING  OF  COHESIONLESS  SOILS 


The  resilient  modulus  of  cohesionless  soils  Is  dependent 

upon  the  magnitude  of  the  confining  pressure  o ^  and  is  nearly  indepen¬ 
dent  of  the  magnitude  of  the  repeated  axial  stress.  Therefore,  it  is 
necessary  to  test  cohesionless  materials  over  a  range  of  confining  and 
axial  stresses.  (The  confining  pressure  is  equal  to  the  chamber  pressure 
less  the  back  pressure  for  saturated  specimens.)  The  following  proce¬ 
dures  should  be  used  fu.  this  type  of  test: 

a.  Use  confining  pressures  of  5,  10,  15»  and  20  psi.  At  each 
confining  pressure,  test  at  five  values  of  the  principal 
stress  difference  corresponding  to  multiples  (l,  2,  3,  5)  of 
cell  pressure. 

b.  Before  beginning  to  record  deformations,  apply  a  series  of 
conditioning  stresses  to  the  material  to  eliminate  initial 
loading  effects.  The  greatest  amount  of  volume  change  occurs 
during  the  application  of  the  conditioning  stresses.  Simula¬ 
tion  of  field  conditions  suggests  that  drainage  of  saturated 
specimens  should  be  permitted  during  the  application  of  these 
loads  but  that  the  test  loading  (beginning  in  Step  f  below) 
should  be  conducted  in  an  undrained  state. 

£.  Set  the  axial  load  generator  to  apply  a  deviator  stress  of 

10  psi  (i.e.,  a  stress  ratio  equal  to  3).  Activate  the  load 

generator  and  apply  200  repetitions  of  this  load.  Stop  the 
loading. 

d_.  Set  the  axial  load  generator  to  apply  a  deviator  stress  of 

20  psi  (i.e.,  a  stress  ratio  equal  to  3).  Activate  the  load 

generator  and  apply  200  repetitions  of  this  load.  Stop  the 
loading. 

£.  Repeat  as  in  Step  d  above  maintaining  a  stress  ratio  equal  to 
6  and  using  the  following  order  and  magnitude  of  confining 
pressures:  10,  20,  10,  5,  3,  and  1  psi. 

£.  Begin  the  record  test  using  a  confining  pressure  of  1  psi  and 
an  equal  value  of  deviator  stress.  Record  the  resilient 
deformation  after  200  repetitions.  Increase  the  deviator 
stress  to  twice  the  confining  pressure  and  record  the  resilient 
deformation  after  200  repetitions.  Repeat  until  a  deviator 
stress  of  5  times  the  confining  pressure  is  reached  (stress 
ratio  of  6) . 

£.  Repeat  as  in  Step  £  above  for  each  value  of  confining 
pressure. 

h_.  When  the  test  is  completed,  decrease  the  back  pressure  to  zero, 
reduce  the  chamber  pressure  to  zero,  and  dismantle  the  cell. 
Remove  the  LVDT  clamps,  etc.  Remove  the  soil  specimen,  and 
use  the  entire  amount  of  soil  to  determine  the  moisture  content. 
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APPENDIX  B:  CHEVDEF  PROGRAM 


INTRODUCTION 

The  CHEVDEF  program  takes  measured  deflections  from  a  deflection 
basin  with  initial  estimates  and  ranges  of  layer  modulus  and  computes 
the  modulus  values  that  best  describe  the  input  deflection  basin.  A 
linearly  layered  elastic  computer  program  originally  developed  by 
Chevron  Oil  Company  is  used  as  a  subroutine  to  calculate  the  stress, 
strains,  and  deflections. 

The  information  provided  herein  is  as  follows: 
si.  Flowchart . 
b .  Input  guide . 
c_.  Example  input. 

d.  Example  output. 

e.  Program  listing. 

FLOWCHART 

A  flowchart  describing  the  logic  of  the  program  is  presented  on 
the  following  page. 


INPUT  GUIDE 

PROGRAM  CHEVDEF 

(Matching  of  Pavement  Deflections 
Using  CHEVRON  Layered  Elastic 
Computer  Program) 

Line  1 

LINENU  I  NPROB  j 


NPROB  =  Number  of  data  sets 

Line  2 


Title  72  characters 

Note:  Line  2  through  Line  13  are  repeated  for  each  data  set. 

Line  3 


LINENU 

ND 

RRD(l) 

1 

1 

1 

1 

RRD(ND) 

ND  =  Number  of  deflection  reasings  (maxinnmi  of  b) 
RRD(i)  i  =  1,  ND  ~  Measured  deflections  in  mils 


Line  L 


LINENU 

NL 

TOI. 

MAXIT 

NL  =  Number  of  variable  layers  for  which  the  modulus  is  to  be 
determined  (not  to  exceed  the  number  of  deflections) 

TOL  =  Tolerance  in  percent  for  stopping  programs  (usually  =  10) 

MAXIT  =  Maximum  number  of  iterations  (usually  =  3) 

Line  5-1  through  5-NL  (one  line  for  each  unknown  modulus) 


LINENU 

ILV(i) 

EMIN(i) 

_ 

EMAX(i) 

ILV(i)  =  System  layer  number  for  unknown  modulus  value  i 
EMIN(i)  =  Minimum  allowable  modulus  for  unknown  modulus  i 
EKAX(i)  =  Maximum  allowable  modulus  for  unknown  modulus 

Line  6 

Title  at  start  of  CHEVRON  data  72  characters 


Line  7 


LINENU 

WGT 

HSU 

NOUTP 

NPUN 

WGT  =  Total  load  applied  to  pavement 
PSI  =  Contact  pressure  of  load 
NOUTP  —  Control  output  :  set  -  0 
NPUN  =  Control  output  :  net  ■“  1 

Line  8  (continue  to  additional  lines  as  necessary  for  all  E's  and  V's) 


LINENU 


NS 

E(i) 

V(i) 

— 

— 

E(NS) 

V(  NS ) 


NS  =  Total  number  of  layers  in  system 
E(i)  =  Starting  modulus  for  layers  i  =  1,  NS 
V(i)  =  Poisson’s  ratio  for  layers  i  =1,  NS 


Line  9 


LINENU 

HH(1) 

HH(P) 

-  •  - 

1 

1 

HH(NS-l) 

0 

HH(i)  =  Thickness  of  layer  i  -•  1,  NS-1 

Line  10 


LINENU 

iR 

**’ .  ** 

RR(1) 

RR(  2 ) 

— 

RR(iR) 

n 

ifi  =  Number  of  radial  offsets  (set  =  ND) 

RR(i)  i  =  1,  iR  =  Distance  to  deflection  readings 


Line  11 


LINENU 


iZ 


ZZ(l) 


iZ  =  Number  of  depth  set  -  1 


ZZ(l)  =  Depth  of  deflection  set  =  0 

Note:  After  determination  of  final  modulus  subgrade,  CHEVRON  is 
called  for  computation  of  stress,  strains,  and  deflection 
at  selected  points.  These  soils  should  be  the  center  of 
granular  layers  and  the  top  of  the  subgrade  material. 


LINENU 

iR 

RR(1) 

RR(2) 

_ 

RR(iR) 


I 

iR  =  tlumber  of  offsets  (f  iR  =  0  returns  for  new  data  set) 
RR(i)  i  =  1,  iR  -  Distance  to  selected  points 


Line  13 


LINENU 

iZ 

ZZ(1) 

ZZ(2) 

RR(iZ) 


iZ  =  Number  of  depth 

ZZ(i)  i  *  1,  iZ  =  Depth  to  selected  points 

Note:  Run  is  now  terminated  or  returned  for  new  data. 

For  devices  with  two  loaded  areas,  use  total  force 
on  one  area  and  compute  radial  distances  for  Line  10. 

EXAMPLE  INPUT 

An  example  input  for  the  PTRF  1A  item  is  listed  below.  Control 
cards  preceeding  input  are  for  the  Honeywell  G635  Computer  in  the  re¬ 
mote  batch  mode  (CARDIN). 


01  0««N  <i  ) 

02 Oil  I  DENT ! R 0SF3 0 0 -  CURT I S 

03  04 ! DPT  I DN ! FORTRRN 

04 05! USE! . GTLIT 

05  04 ! FDRTY ! XREF 

06  0* ! SELECTh! R  0SF3  0 O'CHEVDEF 

0?  OH  EXECUTE 

08  0* !  L I M I TS ! 4  0 » 3  OK  *  *  6K 

090  010  1 

100  Q20PTRF  lft  RR2 0 03  DEFL.  <7000  LB> 

110  030  4*5.903*4. 433?  3. 067 *  1 . 933 
120  040  3  *  10*3 
130  050  1*20000* 700000 
140  060  2*20000* 100000 
150  070  3* 10000» 30000 
160  080  CHEVRON 

170  090  70  00  *  27. 5  08  * 0  * 1  ^  ^  .  , 

180  100  4  *  5  0  0  0  0  0 • 0 . 35  * 5  0  0 0  0 >  0 . 35  * 2  0 U 0 0  >  0 . 4  >  1 0  0  0 0  0  0  *  0 . 

190  110  7. 5* 20. 5? 212 

200  120  4* 4. 5* 12? 24? 36 

210  130  1*0 

220  140  1*0 

230  150  2*7.5*28 

24  OtENI'JDB 
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IXAXPLE  OUTPUT 


An  example  output  for  the  PTRF  1A  item  is  provided  on  the  fol¬ 
lowing  pages.  The  number  of  problems  to  be  solved  is  1  . 


'-••s 0  PTPF  lft  RR2008  DEFL.  •7000  LB.  ■ 


NUMBER  OF  VRPIftBLE  LftYEPS  ft MB  TftPGET  DEFLECT IONS  =  3 


DEFLECTION  PEftblNGS  IN  MIL .3 
:'OS  IT  ION  NO:  1  £ 

DEFLECTIONS:  K..908  4  4-5 

WEIGHTING  FRC TOR:  0.169  o' 28, 


3 

3.  067 
0.326 


4 

1 . 93  3 
0.  517 


VflPIftBLE 
LftYER  NO 
1 

2 

3 


SYSTEM 
LftYER  NO 
1 

2 

3 


VftLUE  OF 
MfiXMUM  MODULUS 
7000 00. 0 
100000.  0 
30000. 0 


VftLUE  OF 

MINIMUM  MODUL  US 

20000. 0 
20000. 0 
10000. 0 


♦♦♦♦♦080  STftRT  DftTfi  FOP  CHEVRON 


THE  PROBLEM  PftRftMETERS  ftPE 


TOTfiL  LOftD. . 
TIRE  PRESSURE.. 
LOftD  PftDIUS. . 

L ft YEP  NO. 

1 

2 

4 


7  0  0  0 .  0  0 
27.51 
9.  00 

MODULUS 

500000. 

50000. 

20000. 

1  00  0  0  0  0 . 


LBS 
PS  I 
IN. 

POISSONS  RftTID 
0 .  35  0 
0 . 35  0 
0. 400 
0 . 5  0  0 


THICKNESS 
7.50 
20.50 
212. 00 

SEMI-INFINITE 


PO.-ITION 

1 


4 


DEFLECTION 
5.5500 
4.2409 
2.3706 
2. 0131 


MEftSUPEB 

5 . 9  03  0 
4.4330 
3.  0670 
1 . 9330 

UM: 
SUM: 


ftBSOLUTE 

ftPITHMETIC 


DIFFERENCE 

0. 3580 
0. 1921 
0. 1964 
-  0 .  08  0 1 


*.  DIFF. 

6.  1 

4.3 

6.4 
-4.  1 


0.8265  20.9390 

12.6527 
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DATA  FOR  DEVELOPING  EQUATION  I  POP  ITERATION: 


1 


NO. 


LAYER 

INITIAL 

CHANGED 

OFFSET 

df flection: 

r*D. 

•5 1 

MODULUS 

MDDULU : 

DISC. 

INI T IAL 

C HANGE  D 

READ  IN 

i 

S  ft  0  0 ft  ft . 

£  0  0  0  0 . 

4.50 

S  .  SS  ft 

14.501 

5.  -?08 

IS.  00 

4.  £41 

5.428 

4.433 

£4.  00 

£.871 

3.  023 

06  r 

36 .  0(1 

£.  013 

2.  102 

1.833 

£ 

50000. 

£  0  0  0  0 . 

4.50 

5 . 55  0 

7.  686 

5. 808 

IS.  00 

4.241 

5.  881 

4.433 

£4.  00 

2. 871 

3.808 

? .  067 

36.  00 

2.  013 

2.331 

1.833 

“j 

J 

£  0  0  0  0 . 

1  0  0  0  0 . 

4.50 

5. 550 

7.507 

5 . 8  08 

IS.  00 

4.241 

6.  037 

4.433 

24.  00 

2.871 

4.608 

3.  067 

36.  00 

2.  013 

3.542 

1 . 333 

PREDICTED  E 

DISREGARDING 

BOUNDARY  CONDITIONS 

5  0  0  06  0 . 

40304. 

21130. 

POSITION 

DEFLECTION 

MEASURED 

DIFFERENCE 

y.  DIFF 

1 

5.8718 

5.3080 

0. 0362 

0.  6 

2 

4.4734 

4.4330 

-0.  0464 

-1.0 

2.3416 

3. 0670 

0.  1254 

4.  1 

4 

1 . 3882 

1.3330 

-0.  0562 

-2.9 

AB 

SOLUTE  SUMS 

0.2642 

8.6555 

ARITHMETIC  SUMS 

0. 7480 

AVERAGES 

0.  0661 

2. 1633 

THE  PI NHL  MODULUS  VALUES  ARE 


500060. 

DEFLECTIONS 

40304.  21130. 

APE  IN  TOLERANCE 

1 000000 

R  2 

VERTICAL 

TANGENTIAL 

RADIAL 

SHEAR 

BULK 

0.  -7.5 

STPE 

-9.397E  00 

4.723E  01 

4.723E  01 

0. 

8.507E  01 

STRA 

DSPL 

-8. 431E-05 

5 . 3  09E  -  03 

6. 797E-05 

6. 797E-05 

0. 

5. 103E-05 

0.  7.5 

STRE 

-9.397E  00 

-8. 452E-01 

-8. 452E-01 

0. 

-1.1 09E  0 1 

STRA 

DSPL 

-2. 135E-04 
5. 909E-03 

6. 797E-05 

6. 737E-05 

0. 

-8. 253E-05 

0.  -28. 0 

STRE 

-2. 01 9E  00 

1.313E  00 

1.313E  00 

0. 

6. 064E-01 

STRA 

DSPL 

-7. 291E-05 
3. 439E-03 

3.  87 IE- 05 

3. 37 IE- 05 

0. 

4.514E-06 

0 .  28 .  0 

STRE 

-2.013E  00 

1 . 701E-02 

1 . 70  IE- 02 

0. 

-1.885E  00 

STRA 

DSPL 

-3. 621E-05 
3. 439E-03 

3. 871 E- 05 

3. 87  IE- 05 

0. 

-1.378F-05 

END  OF  PROGRAM  ♦♦♦♦♦ 
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0  0?  1  0  •'  H8  *  4  ■  ♦  hZ  •  38*  *  »  R8  •  3’8*.  ?  *  »  B8  •  3'-*t«  ?  •  •  f  8  •  ~:  'Mb  • ?  •  • 

o 0*80  0 ' 38*  •  ?  ■  •  R  j  •  3 x*'  ■  •  P  J  1  •  38*  ■  «  F  JO • ?8t  1  *  TITLE'  8 0  ■  « 

o 0?  30  T  E  .  T  ■  n  .  .  FZ  •  1  0 0  •  .  ;  • .  ,  4 «  4  •  ■  .*  C  «  4  •  .  FM  ■  8*  8  '  • 

0  »'?4  O  :*  F'M  •  4  •  4 »  4  •  *  DUMMY • F 4 ' 


H-B 


THi:  PP0GP8M  Comtrih:  hM  I TEPRT  I VE  PROCEDURE  td  C hLC ULhTE 

modulu:  vrlue:  fop  up  td  4  frvement  lh.ep:  ffdm  the  deflec  t  i dm 
eh : 1 m  mercupement: .  the  fpdgprm  u:ec  the  chevpdm  m  -  lryep 
elm: tic  :v stems  fpqgprm  r:  r  :ubpdutime  td  crlculrte  deflect  tom: 
FROM  IMITIRl  RMIi  ITEPRTIVE  MODULI'*"  VRl.UES . 


PPDGPRM  NRME •  CHEVDEF 

CODED  PV:  DP.  MRLTEP  P  BRP*- EP 

POINT  OF  CONTRCT: 

RL DEPT  J.  BUSH.  Ill 
GEOTEC HM I CRL  LHBOPRTOPV 

mrtepmrv:  experiment  strtion 

VICKCPUPG.  M I S  C I S  S I PP I  38180 

COMPUTER:  URTERMRY '  EXPERIMENT  ‘TRTION.  GE*3? 

LHNGUhGE :  FOPTPRN  IV 

DhTE  COMPLETED:  .JULY  1880 

TPECIRL  PEOUIPEMENT:  CRRDIN5  REMOTE  BhTCH  FpQCESCING 
I TOPRGE :  DICC 


I 


M  n*;.  1  0 
0063  0 

0  0  6  r‘0 

0  064  0 
0  065  0 
00660 
00670 
OUt'&O 
00690 
00700 
0  fi  no 
u  0780 
007  30 
00740 
00750 
00760 
00770 
00780 
00780 
00800 
0  08 1  0 
o  08  c'  0 
0  iix  30 
00840 
0  085  0 
0  086  0 
00870 
00880 
00880 
00800 
0081  0 
00880 
00880 
OU840 

0  085  0 
0086  0 
00870 
00880  ' 
00880 
0 1  0  0  0 
0 1010 
01  080 
01  0  30 
0 1  04  0 
01  050 
0 1  06  0 
01070  ] 
01080 
0 1  080 
01  1  00 
0111  0 
01  180 
0  ‘  '  3  0 
"1  140 


DATA  i:  QNF'AF  N 

READ -5. £10-  LI MENU.  NF'FOB 

I.IFITF  '6*1  05  '  NF  FOB 

DD  8888  NPR= 1 • KFFDB 

MR  I  TE  -6.  1818  ■  JIFF 

CONTINUE 

READ -5. 100m- T  EXT  1 
MP I TE -6*  £0  0 • TEXT  1 
F’E AD  -  5  •  £  1  0  -  L1NEMU.ND.  -PP 


L I MENU* NL  * TOL.MAXIT 


.  EMAX  <  I  ' 


PEHP-5. 1 0  0  o - T  EXT  1 
MR  I  TE  -6*800  •  TEXT  1 

READ  -5.810-  L  1  HE  MU.  HP.  -FPIl.  1  -  .  1  =  1 .  NP  ■ 

:um  =  0. 

i;OUMTF=l 
1  CONTINUE 
I  TEF  =  o 

FEfiD-5.810-  LINE  NU . NL  > TOL.MAXIT 

po  5  I  =i.np 

l-i  ■  I  -  =  1 .  FPP  -  I  - 
5  CONTINUE 
DO  8  I  =  1  <  NL 

PE  fiV  -5.810-  LI  NENU .  Nl-l .  Eh  I N  <  I  >  .  EMAX  <  I  ' 

8  ILV ' I • =NW 

MP I TE - 6  *  6  0  0  -  NL 

MR I TE  <6 -618 '  -1.1  =  1. NP> 

MP  I  TE  -  6 . 68  0  :•  <  FPP  -  I  -  .  1  =  1 .  NO) 

MP I TE -  6 . 68 1 >  CM <|) . 1  =  1. HP  - 
l-lPITE  -6.  780  - 
WRITE  <6.650) 

PD  58  MS= 1 . NL 

MP I TE  <6 » 66  0>  MS  *  I L V  CMS) . EMAX -MS > . EMIN  CMS) 

8  CONTINUE 
CALL  CHEVRON 
WR  I  TE  -  6 . 83  0  - 
CONTIN  =  .TRUE. 

Z  UM  =  0. 

•  IJMP  =  0 .  0 

a:um=o. o 

PD  8  I  =  l . HP 
err  =  PPD -  I • -DEE  <I) 

PERC=  -ERR  - RPP < I • ) 41 00 
IUM  =  SUM  ♦  MBS -ERR) 

:UMP  =  SUMP  ♦  MBS  <PERC - 
ASUM  =  MSUM+PEPC 

WRITE -6.880- I . PEF -I) . RPB 1  I ) .ERR.PEPC 
PF ' 1 . 3  -  =  PEF  -  I ) 

CONTINUE 

IF  <  SUMP. AT. TDL'  CONTIN  =  .FALSE. 

WRITE  -6.840-  SUM. SUMP. ASUM 
IF  -.CONTIN-  GOTO  88 
C OUNTR=  3 
I  TEP= I TEP  + 1 
PO  15  I  =  1 . NL 
^  =  ILV- I ' 

ELCI.3*  =  ALOG 1 0  < E - K ' > 

CONTINUE 

MR I TE < 6 . 760 )  ITER 
PO  60  JX* 1 . NL 
K  *  I L  V  -  JX) 

IF-ITEP.EO. 1 ) GO  TO  £1 
I F  <C  -  JX) . EO . C TEMF C JX-  - 

S  C  TEMP  -  JX>  «< ftLOS  10  -EMAX  CJX>  '  -t-ALOGl  0  <EMIN  <  JX' 
EL  <  JX.  1  '  =  <  3 ♦»:  -  JX  >  *C  TEMP  >  JX  -  >  .-  4 


m  isn 

IF  ‘ftB  '  E  •  )  ■  -EMfiX  ■  j: :  >  •  .  L T.  l oo. DR. hb .  ■  E  • )  ■  -FH1N  •  i: . 

■  '.LI.  100. GO  T  0 

"  1 1  -■  i'i 

I F  ■  i;  ■  J. :  •  .  L  T .  t  TEMP  •  j: :  •  •  GO  TO  23 

"I  l  ,'n 

i  f  ■  m  ‘  i.Ki .  l  t  .  o  ■  fl  •  j: : •  i  •  =  •  3*c  ■  .ix>  +mlci&  i  o  ■  EMftx  •  j: :  •  •  • 

4 

111  lx  II 

GO  TO  24 

Ml  190 

83 

IF  •  H:iJM.i3T.  0>EL  ■  JX.  1  '  =  •  3*C  '  JX>  ♦ftLQG  1  0  •  EM I N  •  JX 1  '  > 

■  4 

N  1  M  0 

24 

E  '  ►  •  =  1  '.'♦♦EL  JX.  1  ' 

0 1 2  1  0 

GO  TO  40 

diced 

21 

IF  •  El-IJM.  LT,  0  '  GO  TO  30 

012  30 

E  <►  >  *EHIN  <  J>: » 

01240 

GO  TO  40 

012®.  0 

30 

E  <Y  ■ =EMfiX ■ JX  * 

0 1  £>■  0 

40 

EL  ■  JX.2> *fiLDG10'E  ‘>  )> 

01270 

CRLL  CHEVRON 

01280 

DO  50  M*  =  1 « NF 

01280 

DF  <  M*  »  2 •  =  DEF <YY> 

0 1  3  0  0 

45 

CONTINUE 

0 1310 

*  C.K  ►  *  JX  >  —  (  D F  '  ►  K  *  2)  —  DF  i  K  K  9  3>  )  /  *.EL  '•  JX  •  d  *  ~EL  •  JX  •  3.*  • 

01  320 

M  <Y  Y .  JX  >  =DF  <YY  .  2>  -EL  <  JX*  2>  *S  <YY  »  JX) 

01  530 

50 

CONTINUE 

01  340 

ET3  =  1 0**EL  <  JX*  3) 

01  350 

WRITE  <6. 770) ILV<JX> *ET3*E <ILV < JX) > »fiXl < 1 > . DF < 1 . 3 

•  .DEF' 1  * .RRD'  1  ' 

01  36  0 

IF'ND.EO. 1)  GO  TO  39 

01  370 

DO  34  JS  =2*  ND 

01380 

WRITE  <6.  780) fiX  1  <  JS>  *  I'F  <  JS*  3)  *  DEF  < JS)  *  RRD <  J$> 

01  380 

34 

CONTINUE 

0 14  00 

39 

CONTINUE 

0141  0 

WRITE '6. 810) 

01420 

E  '>  >  =  1 0**EL  •  JX*  3 » 

014  30 

60 

CONTINUE 

01  440 

DO  65  M<  =  1 • ND 

01450 

65 

CONTINUE 

0 1 46  0 

DO  69  Y  Y  *  1  *  ND 

01470 

88 

CONTINUE 

01480 

IF  <NL.NE.ND>  GO  TO  101 

01480 

DO  120  I  *  1  * NL 

01500 

C  a  >  =  RRD  •  I  >  -ft  <  I .  NL  * -i  'I*  NL>  *EL  « NL .  3- 

01510 

DO  110  J  =  1  *  NL 

01520 

C  •  I  >=C  n  <1*  J*  «EL  <  J,  3, 

015  30 

B'l.Ji  =  3<I*J> 

01540 

BTEHF  <  I  »  J  :■  =  B  •  I  .  J> 

01550 

110 

CONTINUE 

0 1 56  0 

CTEMP • I>  =  C • I» 

01570 

120 

CONTINUE 

01580 

GO  TO  79 

01580 

101 

CONTINUE 

0 1 6  0  0 

DO  30  I  =  1 « NL 

0161  0 

c  <  I  •  =  0.0 

01620 

DO  70  J  =  1  *  ND 

0 18  3  0 

C  •  I  >  =  C  'll  *S  '  J.  I  '  ♦  <PPD  •  J*  -ft  i  J.NL  '  '  J.  NL  >  *EL  •  NL  • 

3  •  •  *u  •  j » 

01640 

DO  67  j:  *  1 , NL 

01650 

C  'I  i  =C  'I '  ♦  I  ‘  J.  JS  •  *S  '  J.  I  '  *EL  '  JS »  3  >  *W  •  J  • 

01660 

6? 

CONTINUE 

01670 

70 

CONTINUE 

01680 

DO  68  J3  =  1  *  NL 

01680 

B ■  I  *  JS »  =•  0. 

01700 

DO  71  J  *  1 « ND 

01710 

B ■ I «  JS i  *  B'l.Ji'  ♦  S' J* I**S • J* JS'*W' J' 

01720 

71 

CONTINUE 

017  30 

BTEHP'l* JS*  «  B'l.Ji* 

01740 

68 

CONTINUE 

fl 


B-10 


n.»  ru 


i_  •  i 


ul  T50  L  TEMP '  I  • 

nl  760  :?0  CONTINUE 

'  .1770  7  ?  i;  OUT  I  HUE 

..  17>'0  DO  78  MM=  1  •  ML 

ii  1  780  78  C  OUT  I  HUE 

01800  OEtLL  I  IMEC"  £•  C  •  ML  •  k  EP»  4  • 

i‘  181**  UFITE '6.  :: 0  1  * 

01  88 1 1  501  FOPMPT  •  tH  • 

0 1  8  5  0  ME  I  TE  1  0.  •  ?■  ft  1  • 

0 1  84  0  UP  I  TE  ‘  6 «  5  08  ' 

01850  508  FOPMhT  .  1H  .  PEE II IC  TED  E  Ii I  CPEGmPD I Mb  POUNDPPY  CONDITION!  • 

0 1 8*  o  £'0  7  5  j; :  =  t .  ML 

0 1 87  0  hTEMp  .  j: :  •  =  l  o**r  ■  j:  ■ 

01880  75  c OH T I HUE 

01880  ME  I  TE  *6.8-51  '  **TEHP*  .111  -  .  JX»1»ML> 

01800  IF"EP.ME.  0"30  TO  80 

01810  EECflLC  = .  FhL  5E. 

01880  DO  84  I=1.HL 

018  50  I F  ■  PEC FiLC  '  '50  TO  84 

0 1  84  0  pmp:  :  =  HLOO 1  0  ■  EMPM  •  I  >  > 

0 1  85 0  HM I M  =  ML 013 1 0  *  EMIN ■ I ■  ■ 

01880  I F • C  *  I • . 6E. 0001- . AND. C • I>.LE.  <AMAX4-.  0001 ■  "3D  TO  84 

01870  FECmLC - . TEUE. 

0188  0  I'D  88  k  = 1 <  NL 

o  1 88  0  IF"  .  EG' .  I  1  bO  T 0  I;. 6 

08  0  0  0  C  "  ■  =C  TEMP  •>:  '• 

0801  0  86  IiO  81  L  =  1  <  ML 

080  £  '  L » ►  7'  =  BTEME  •' L « k > 

030  81  C DM T I HUE 

0804O  88  CONTINUE 

08  05  0  IF  "!  •:  I>  .LT.RMI N> C  < I > =RM I N 

08060  IF  "5  ■  I  •  .  i3T.  Rflfi:*:  ■  C  <  I  >  =FtMElX 

08  0  7  0  C  TEMP  <  I  :>  =C  '  I  > 

08080  DO  83  1=1' ML 

08080  B 1  I » !■ =  0 . 0 

08100  PTEMP ' I >  k  '  =E  •  I  ■ k> 

081  10  83  CONTINUE 

08180  B' I' I  1=1.0 

081  30  BTEMP | I • I i =  B • I ' I > 

08140  84  CQMTIMUE 

08  150  I F  •  EEC  FiLC  .*  >30  T □  78 

08160  BO  85  1=1 'ML 

08170  J=IL  V i I ' 

08180  TEHE=tO**C <I> 

08  !  8  0  TEMP  1  =  1 0**EL • I «  3 ■ 

08800  ETOL -  I  - =HB! -  •  'TEMP-TEMPI •  'TEMPI >  *1 00) 

08810  E  •  J  • =TEMP 

0888  0  CTEMP • I ' =EL • I  *  3 • 

08830  85  CONTINUE 

08840  ChLL  CHEVPON 

08850  I  l.iM  =  0. 

08860  !UMP  =  0,0 

08870  Ft'l.lM  =  0.0 

08880  C DMT  IN* . TPUE . 

08880  i.tPITE  '6'830‘ 

08300  BO  88  I  =  1*ND 

08  310  EPF  =  FED < I • -BEF  > I  • 

08 380  FEFl  =  •  EPP  PPI1  ■  1  •  •  •  1  00 

■8  3  30  :ur  =  :um  ♦  mb : .  epf  ■ 

08  34  0  ;unp  =  :ump  ♦  hb : ■ fepc  < 


B-]  1 


I.»C  .-I  'll 
f»c*  ii 


05  ?  r*  f i 

05  4  0  »i 

0£  4  1  fi 

03450 

054  3  0 

0544  0 

05450 

0546  0 

03470 

05430 

05490 

05500 

055 1  0 

05550 

055  3  0 

0554  0 

05550 

05530 

05570 

05530 

05590 

05300 

053 1  0 

05350 

053  30 

0534  0 

05350 

0533  0 

03370 

0533  0 

05390 

037  00 

057 1  0 

05750 

05730 

05740 

05750 

0573  0 

05770 

05780 

05790 

05800 

058 1  0 

05850 

058  30 

05340 

05350 

05830 

05870 

05880 

05390 

05900 

059 1  0 

05950 

05930 

0594  0 


Up  I TE  *3.850  *  I . OEF  *  I > . PPf *  I  *  « EPP.PEPC 
DF  •  I  . 3 ■ =  DEF *  I  * 
h:um  =•  h:ijm  ♦  ferc 
88  COtlT  I  HUE 

IF  • SUMP.GT. TOL ■  contin=.fhl :e. 
per:um=:um  nd 
PERSUMF=SUMP  ND 
UP  I TE  ' 3 • 84 0 >  SUM, SUMP. AS UM 
UP  I TE ' 3 . 7 1 5  <  PEPSUM . PEP! UMP 
PEITEP  =  . TPUE . 

DO  87  1  =  1 «  NLAYS 

IF 'ETOL ' I >  . 6T. TOL>  PE  I TER=. FALSE. 

87  CONTINUE 

IF  'CONTIfO  GOTO  89 
IF 'PEITEP *60  TO  39 
IF  *  I TEP . LT . MAX IT) GO  TO  3 
WRITE <6, 9000) >E'I)  * I*1»NLAVS> 

WPITE <3. 31 0> 

UP  I  TE  *'3.  303* 

303  FOPMRT'IH  .REACHED  MAX  NO  DF  ITERATIONS') 

GOTO  10 

39  WRITE <3,9000> *E  '  I) , 1  =  1 , NLAYS) 

I F  *  C  ONT  IN  •  I..IP  I  TE  <3,  3 04 > 

304  FOPMRT'IH  ,  DEFLECTIONS  RPE  IN  TDLERRNCE ' .) 

IF ( PE  I  TER )  WPITE <3*3  05  > 

305  FOPMRT  < 1 H  ,  '  CHANGE  IN  MODULUS  VALUES  RPE  IN  TOLERANCE) 

10  C0UNTP=99 

CALL  CHEVRON 

GO  TO  9999 
90  UP  I  TE  *'3.87  0'* 

9999  CONTINUE 

WPITE *  3.850* 

ST  DP 

105  FOPMRT'IH 1.1  OX,  THE  NUMBER  OF  PROBLEMS  TD  BE  SOLVED  IS'. 

4  8"  ’ ,  I  3 *  *■  ■  ■  * 

380  FORMAT '5X» 15. 3F 1 2 . 4 ,  F 1  0 . 1) 

830  FORMAT*  ,5X,  POSITION'  .  8X.  DEFLECTION'  ,  3X,  'MEASURED' , 3X, 

4-  DIFFERENCE'  ,  3X»  X  IHFF.  > 

9000  FOPMRT  3X»  THE  FINAL  MODULUS  VALUES  APE',  .'V, 

4  4  *  F 1 0. 0, 4X> ) 

500  FORMAT <1H  » A80 > 

810  FORMAT 'V* 

840  FORMAT*-,-.  T30,  IN  THE  A  MR  TP  I X  .  '  ,  T  1  3 «  1  '  ,  T  3  0 .  8  '  •  T  48 ,  3  . 

4T38.  '  4  • 

55 1  FORMAT  5X ,  4  •  F I  0 .  0 ,  3X*  * 

85 0  FORMAT  '8X ,  4  CE 1 3 ■  3 , 3X)  .:* 

830  FORMAT*:  v,T30*  IN  THE  S  MATRIX  .  *v  ,  T13,  '  1  T30,  T48,  '  3  •» 

4  T 38 ,  4  '  * 

870  FORMAT'  THIS  MATRIX  HA:  NO  SOLUTION) 

290  FORMAT  '  '  T30«  '  IN  THE  B  MATRIX '  .  .  Tl  3,  '  1  ,  T30,  '5'  ,  T48.  '3 ' , 

4  T 38 «  4  '  « * 

300  FORMAT  ‘".T30.  IN  THE  C  MATRIX  .  ' ,  Tl 3,  1  ,  T30,  5'  , T48,  3  » T38.  4 

!*':  4  •  • 

310  FORMAT *T85.  THE  MODULUS  VALUES  APE  NOT  WITHIN  TOLERANCE') 

1000  FORMAT 'A30> 

300  FORMAT '5X.  NUMBER  OF  VARIABLE  LAYERS  AND  TARGET  DEFLECTIONS  = 

"-I2,'  * 

319  FORMAT < IX. //. 15X.  DEFLECTION  READINGS  IN  MILS','. 

4  POSITION  NO:  , 3X , 4 • 9X , 1 1 , 5X * ) 

32 0  FORMAT '  DEFLEC  T IONS : ' . 3X.  4 ■ F 1 2 .  3 ■ ) 


02950  62 1  FORMAT  <  WEIGHTING  FACTOR:  *4*F12.  3>  > 

02*60  650  FORMAT-  ,5X,  VARIABLE  SYSTEM  •  9X*  VALUE  0F*12X.  VALUE  OF '  »  '  * 

02970  ?-5X,  LAYER  NO' *  3X»  LAYER  NO ' * 5X »  HAXNUM  MODULUS  *5X»  MINIMUM  MOSUL 

02980  % U 3 '  • 

0299 0  66 0  FORMAT <9X* 12,  9X ,  1 2 ,  1  OX ,  F 1  0 .  1  «  1  OX , F  J  0 .  1  > 

03000  710  FORMAT <5 </> «20X,  FOR  EQUATION  A  ♦  J:  ♦  E  =  DEFLECTION' » /* > 
03010  720  FORMAT  <.5 *v>»23X»  FOR  THE  EQUATION  [PI  ♦  CE]  =  £CJ',  -  > 

03020  ’61  FORMAT  *80  <'!«♦>  ,5X,  DATA  FOP  DEVELOPING  EQUATIONS  FOR  ITERATIONS'. 

03030  %'  NO.  «2X*I2> 

03040  760  FORMAT <1H1 .5- •> *5X«  DATA  FOR  DEVELOPING  EQUATIONS  FOR  ITERATION" 

03050  *  NO.  • I2«  ss/HXt  LAYER  INITIAL  CHANGED  OFFSET', 10X. 

03060  t  DEFLECTIONS' ,  -'5X, ' NO. ' ,5X,  MODULUS  MODULUS' *5X» '  DISC.  INITIAL 
03070  y-  CHANGED  READINGS' 72 <’'♦•')  > 

03080  770  FORMAT '6X« I2,4X,F8. 0, 3X.F8. 0, 4X,F6. 2, 3X, F6. 3, 6X, F6. 3, 3X.F6. 3> 

0 3 090  73 0  FORMAT  •  T  36 » F6.  2 ,  3X,  F6 .  3, 6X * F6.  3, 3X ,  F6.  3) 

0 3 100  790  FORMAT  «  -  ✓  -  > 

031  10  310  FORMAT -72 • 

03120  715  FORMAT r2SX« "AVERAGE: "»2F1 0.4) 

031  30  840  FORMAT, 23X* "ABSOLUTE  SUMS " *  2F 1 0 . 4 » / » 22X , 

03140  ARITHMETIC  SUMS ” , 1  OX, FI  0. 4> 

03150  850  FORMAT ,  1H  ,26H*****  END  OF  PROGRAM  ♦♦♦♦♦> 

03160  1212  FORMAT <1H1 »///, "PROBLEM  NUMBER  ",  14, s//s> 

03170  END 

03180  SUBROUTINE  CHEVRON 

03190CR040  65P040  N-LAYER  ELASTIC  SYSTEMS  PROGRAM 

032 00C  CALCULATING  STRESSES,  STRAINS,  AND  DEFLECTIONS 

0321 0C191MN  ♦♦♦♦♦♦  MAIN  ROUTINE  -  N-LAYER  ELASTIC  SYSTEM  ♦♦♦♦♦♦ 

0322 OC 


0  32  3  0 

COMMON 

RMCOY,  RP(10>  , 

ZZ<10>, 

E  <5)  , 

V  <5)  « 

HH  <4  >  , 

0324  0 

& 

H  <4>  , 

AZ  <396) , 

A <396, 5) , 

B  6396, 5)  , 

C  <396,5) , 

03250 

& 

D<396,5> , 

AJ  <396) , 

RJ1 <396) , 

R  J  0  <  396)  * 

TITLE <20) , 

0  326  0 

& 

TEST  <1 1)  » 

BZ<100)  , 

X  <5,4,4) , 

SC <4> , 

FM ■ 2 , 2) , 

0  3270 

PM  <4,4,4) 

,R, 

Z, 

AR, 

NS, 

03280 

S' 

N, 

L» 

ITN« 

RSZ, 

PSP, 

03290 

& 

POM, 

RMU, 

SF, 

CSZ, 

CST , 

0  3  3  0  0 

% 

CSR, 

CTR, 

COM, 

CMU, 

PS  I, 

0  3  3 1  0 

S' 

NLINE, 

NOUTP, 

NTEST , 

I, 

I  TN4, 

03320 

S- 

K , 

LC, 

JT, 

TZZ, 

PR, 

03  3  3  0 

S' 

PA« 

P, 

EP, 

TIP, 

TIM- 

0  3  34  0 

& 

Tl, 

T2, 

T3, 

T4, 

T5, 

03  350 

& 

T6, 

T2P, 

T2M, 

WA, 

BJ1, 

03360 

%: 

BJO, 

ZF, 

SZ1, 

SZ2, 

SGI, 

03370 

% 

SG2, 

PH, 

PH2, 

VK2 , 

VKP2, 

03  380 

& 

VI-4. 

VPP4, 

VKK8, 

RDT , 

RDS 

03 39 0  COMMON  OPT.'COUNTP , DEF <  4) , NLAYS 

03400  INTEGER  COUNTR 

03410  DATA  ASTER,  FERD'4H****«  4H....X 

03420  IF  -C0UNTR.EQ.99)  GO  TO  7 

03430  IF (COUNTR. GE.2)  GO  TO  26 

03440  CALL  FLGEOF  <41»NFILE> 

03450C 

0346 OC  ♦♦  COMPUTE  ZEROS  OF  JI <X>  AND  JOOO.  SET  UP  GAUSS  CONSTANTS  ♦♦ 

03470  If  =  ITN+1 

03480  DO  2  1=7, K, 2 

03490  T  =  I '2 

03500  TD  =  4.0*T  -  1.0 

03510  2  BZ f I  *  =  3. 1415927*, T  -  0.25  *  0. 050661- TD 

03520  %  -0.053041  TD**3  ♦  0. 262051 'TD**5> 

03530  DO  3  I=8,ITN>2 

0354  0  T  =  ■  I-2.»  -  2 


B-13 


0  355  0  TP  =  4 .  0*T  ♦  1.0 

035ft0  3  PZ  •  I  '  =  3.  1 4 159£7*  <  T  ♦  0.£5  -  n.  151985  TB 

0  3570  +  0.  0153 '3'?  T  B  •  *  3  -  0.345370  TIi**5> 

i)  if.  ft  0C 

(13590  10  REAP  •  5.  39ft'  TITLE 

03ft 00  39ft  FORMAT ■ ft Oft 4 ■ 

0  3ft  10  IF  'riFILE.EC'.  1  >  6D  TD  9999 

0  3 ft  0  310  FORMAT  '  ft  0A4  ■ 

03  ft- 30  REAP  ■  5  •  399.  L  I  MENU*  WGT  *  F"3  I  *  NOUTP*  NFUN 

0  3ft 4  0  31  1  FORMAT  '£FI£.0*  lift*  113 

03  ft- 50  REAP -5*  399'  L  I  MENU*  NS.  <  E  <  I )  *  V  ■  I  ■  *  I  =  1,  NS) 

0  3  ft- ft-  0  ML  A  Y  3  =  NS 

0 3ft 7 0  301  FORMAT • I 3*F9. 0*9Fft. 0' 

03ft ft 0  399  FORMAT • V • 

0  3 ft- 9 0  1  FORMAT  •  ftX *  1  OFft .  0  • 

03700  1001  CONTINUE 
03710  N  =  u:  -  1 

037c! 0  PEAD>5*  399'  LINENU*  <HH <  I  >  .  1=1 *N> 

0  3 7  3 0  3 1  3  FORMAT  ■  1  OFft .  0 ' 

0  3  74  0  7  PEFiP  <  5  *  399  •  L  I  MENU  *  IP*  <PP  <  I  >  >  I  =  1 »  I P ) 

03750  IF- IP. EC'.  O'  PETUPN 

0 3 7 ft. 0  £001  CONTINUE 

03770  PEAP'5*  399'  LINENU*  12* <22 < I> *  1  =  1 . IZ> 

037ft 0  IF  •  C OUNTP .  EC*.  99*  GO  TO  £ft 

0 3790  3001  CONTINUE 

03800  Ftp  =  3 OPT  'M6T/'  •'  3.  14159#PSI )  > 

03810  NLINE  *  17+NS 

038£0  NPAGE  =  1 

03830  MPITE 'ft. 350)  'TITLE  • I> *  I=1>£0' 

03340  350  FORMAT  ■ 1H1  - 1H0* 5H**»**» £0A4> 

03350  MPITE 'ft* 35 1> 

0  33ft  0  351  FOPMFiT  • 1H0* £3X. £GHTHE  PRDBLEMPAPAMETEPS  APE//) 

03870  MPITE' ft.  35£>  MGT-PSI.flP 

03830  35 £  FORMAT  ■ 1H0*  5X*  12HT0TAL  LOAD..*  SX.  F10.2.  5H  LBS./ 

03890  -i  1 H 0 *  5X *  15HTIPE  PRESSURE. . .  5X*  F10.£*  5H  PSI.- 

03900  i  1H0.5X.  1  3HLOAP  RADIUS.  .  *7X*F10.£*5H  IN./--) 

0 39 1 0  UP  I TE ■ ft  *  353) 

03950  353  FORMAT • 1H0. 5X. 9HLAYEP  NO. * ftX. 7HMQDULUS * ftX* 14HP0ISSDNS  PATIO 

03930  t  ftX »  9HTH I CK NESS  ) 

03940  MPITE ‘ft* 354 ■  <  I  * E <  I  '  *  V ■ I  ' *  HH < I ) *  1*1. N- 

03950  354  FOPMFiT  •  1 H 0 » ftX *13*  9X •  F  1  0 .  0 *  1  1 X * F5 .  3 *  1 4X , Fft- ,  £ > 

039ft 0  MP I TE ■ ft- *  349 •  NS.  E'NS>*  V -H S > 

03970  349  FORMAT  • 1H0. ftX* I  3* 9X* F 1 0. 0* 1 IX* F5. 3* 1 IX. 14HSEMI-INFINITE  //, 

03930  I F  < C OUNTP . NE . 99)  60  TO  £7 

0399  0  MP  I  TE  '  ft-  *  34ft  :• 

04000  348  FORMAT  • 1H  *  IX* 1  HP*  3X* 1H2 » 1 IX* 8HVERT ICAL • £X. 1 0HTAN6ENT I AL * 

0401  o  *  3!!  •  ftHPAP  I  AL .  ftX  *  5H  SHEAR  *  6X  .  4HBULK  > 

04 0£ 0  £7  CONTINUE 

04030C  ♦♦  A P JUST  LAYER  DEPTHS  ♦♦ 

04 04 0  H • 1 ■ =HH ■ 1 . 

04050  PO  £5  I =£. N 

04 Oft 0  £5  H ' I ' =H ■ I - 1 > +HH ' I • 

04070  CALL  CHECK  ■ 1* 

04080  £ft.  CONTINUE 

04090  IF  ■  C OUNTP.  EC'.  99  ■  MPITE  'ft*  348* 

04100  IPT=0 

04 110C  ♦♦  : TAPT  ON  A  NEW  P  ♦* 

04 1  £ 0  100  I  F  T r I p T  + 1 

04130  IF  ■ IPT  IP'  PETUPN 

04 1  4 ft  1  05  P=PP  ■  IPT  ■ 


B-l’i 


t 

i 


i 


04 1  5  0 

DO  31  I  =1,12 

04  1  *  0 

DO  31  J*1.N 

04  170 

T2  =  ABS  <H < J.>  -  22- I>* 

04  ISO 

IF  ■  T2  —  .  i.< 0 0 1  .>  32,  32,31 

04  1 9  0 

32 

1  22  •  I  >  =  -H  ■  J  > 

0420  0 

31 

CONTINUE 

042 1  0 

IF < COUNT*. NE.99>  80  TO  39 

04220 

WRITE <6,  355> 

04230 

NLINE  =  NLINE-t-1 

04240 

355 

FORMAT • 1H  } 

04250 

39 

CONTINUE 

0426 OC 

«« 

CALCULATE  the  partition 

04270 

CALL  PART 

0428 OC 

♦  ♦ 

CALCULATE  THE  COEFFICIENTS  ♦♦ 

04230 

DO  125  1=1,1 TN4 

04  3  0  0 

P=A2 I) 

04310 

1  07 

CONTINUE 

04  320 

IF ■ NS.GT. 5>  GDTO  108 

04330 

CALL  COES  < I ) 

04340 

•30  TO  109 

04350 

108 

CONTINUE 

043S0 

CALL  C015< I> 

04  370 

109 

IF  •*•  115,115,110 

04  39  0 

110 

PR  =  P*R 

04  390 

CALL  BESSEL  -0,PR,Y> 

044  00 

RJ0'I>  =  Y 

044 1  0 

CALL  BESSEL  <1«PR,Y> 

04420 

PJ1 • I>  =  Y 

044  30 

1 15 

PA=P#AP 

0444  0 

CALL  BESSEL  '1,PA,Y> 

0445  0 

A  3  •:  I  >  =Y 

04480 

CALL  CHECK  <2> 

04470 

125 

CONTINUE 

04480 

195 

I2T=0 

0449 OC 

♦♦ 

START  ON  A  NEW  2  ♦* 

04500 

200 

I2T=I2T*l 

04510 

IF • I2T-I2>  205*205,100 

04520 

205 

2=ABS  ' 22  <I2T> ) 

04530 

IF  •  NLINE  -  54  >  207,206,206 

04540 

206 

NPAGE  =  NPAGE  ♦  1 

04550 

NLINE  =  8 

04560 

207 

CONTINUE 

0457 OC 

•  • 

FIND  THE  LAYER  CONTAINING  2  ♦♦ 

04580 

T22  =  0.0 

04590 

DO  210  J 1  =  1 « N 

04600 

J=NS- J1 

046  1  0 

IF • 2-H • J • >  210,215,215 

0462  0 

210 

CONTINUE 

04630 

L  =  1 

0464  0 

GO  TD  34 

04650 

215 

L  =  J*  1 

04660 

IF  ( 22 •  1 2T '  »  33,  34,  34 

04670 

33 

L  =  J 

04680 

T22  =1.0 

04690 

34 

CONTINUE 

04700 

CALL  CALCIN-IRT' 

04710 

IF  1  T Z* 2’  36*  36,  35 

04720 

35 

22-I2T,  =  -22 • I2T> 

04730 

I2T  =  I2T-1 

0474  0 

36- 

CONTINUE 

B-15 


ii4r-.ii  ao  in  cim 

114 7  A II  CALL  EXIT 


i  >  4  770 

04  780 

04  790 

04  ft  0  or 

:  top 

END 

BLOCK  DATA 

048 1  0 

C  QMMON  RMC  0  Y  RR  <10.. 

ZZ  •  1  n .  . 

E  .5"  • 

V  •  5  ■ . 

I 

X 

44 

04820 

H  <4  ■  . 

AZ  <  396 • • 

A ■ 396. 5 > . 

f  *  3'?t.  •  5  * 

>  C  1  396- .  5 >  . 

04830 

£ 

D  1  396 » 5  '•  • 

A J >  396  >  « 

R J1 <  396 » . 

P JO  <396) 

»  TITLE <20< 

04840 

£ 

TEST' 11 *. 

BZ'IOO'. 

X<5»4»4 > . 

SC ■ 4> • 

FM<2»2> . 

04850 

V 

PM'4.4.4> 

,P. 

Z . 

ftp* 

NS. 

04860 

u 

N. 

L. 

ITN. 

PS2. 

RSR. 

04870 

$ 

POM. 

PMU. 

SF. 

C3Z. 

CST  • 

04880 

& 

CSR. 

CTP. 

COM. 

CMU. 

PSI. 

04890 

•/. 

NLINE. 

NOUTP, 

NTEST . 

I. 

I TN4 . 

049  O0 

0 

(•:. 

LC. 

JT, 

TZZ. 

PR. 

0491  0 

PA. 

P. 

EP. 

TIP. 

TIM. 

0492  0 

%' 

Tl» 

T2. 

T3. 

T4. 

T5. 

04930 

& 

T6. 

T2P. 

T2M. 

1.1  A . 

BJ1. 

0494  0 

* 

BJO. 

ZF. 

S2I. 

SZ2. 

36 1 . 

04950 

& 

362. 

PH. 

PH2. 

VK2. 

VKP2 . 

04960 

04970 

•i 

VK4 . 

DIMENSION  ZB<6> 

VKP4 » 

VKK8. 

PDT  • 

RDS 

04980  EQUIVALENCE  < B  Z • ZB  > 

04990  DATA  ZB-  0.  0.  1. 0. 2. 4048. 3. 331 7. 5. 5201 » 7.  0156'' 

05000  DATA  I TN  46/. I TN4  '184/ 

05 0 1 0  END 

05020  SUBROUTINE  BESSEL <NI . X I . Y> 

05030C  ♦♦♦♦♦♦SUBROUTINE  BESSEL  -  N-LAYER  ELASTIC  SYSTEM  ♦♦♦♦♦♦ 

05  04  OC 

05050  DIMENSION  PZ<6> »QZ <6> . PI <6>  .Q 1 <6> .D<20> 

05060  DATA  PZ-'l . 0E0. - 1 . I25E-4. 2. 871 0938E-7. -£. 3449658E-9. 

05070  iL3.9806841E-U.-l. 1536133E-12/.  QZ/-5.  OE-3. 4. 6875E-6. 

05080  *•-£.  3255859E-8.  2. 8307087E-1 0.  -6. 3912096E-12.  2.  3124704E- 18 

05090  S  RI  -'  1.0E0.  1.875E-4.  -3. 6914063E-7.  2. 771 3232E-9. 

05100  S-4.5U4421E-11. 1 . 2750463E-12  Ql/1.5E-2»  -6.5625E-6. 

05110  s  2.8423828E-8.-3.2662024E-10.  7.  143U66E-12.  -2. 5327056E- 1 3-  . 

05180  &  PIS. 1415927' 

051  3  OC 
05 1 40C 

05150  9  N  =  HI 

05160  X  =  XI 

05170  IF  ‘  X-7. 0 >  1 0 » 10* 160 

05 1 8  OC 


05190 

10 

X8=X'£. 0 

05200 

FAC  =-X2*X8 

05210 

IF  <N>  11.11.14 

05220 

1 1 

C  *  1 . 0 

058  30 

Y=C 

05240 

DO  13  1*1.34 

05250 

T*I 

0526  0 

C  =FAC  ♦C  < T ♦T • 

05870 

TEST=ABS  <C>  - 

05880 

IF  •TEST'1  17.17 

05890 

18 

Y*Y+C 

05  3  0  0 

1  3 

CONTINUE 

05  3 1  0 

14 

C  *X2 

05:20 

Y*C 

0*. '  <  n 

DO  16  1*1.34 

0T<  :<40 

T*  I 

B-l  6 


»*i5  JSO 

C  =F HC  *C -  • T* ■ T  +  l . O  '  • 

05  ir.  i  i 

te:t=hb:  .c*  -  io.o** 

n*:.  370 

IF  ■  TE  ;  T '  17-17.15 

380 

15 

v=v«-c 

05  3  '30 

16 

CONTINUE 

05400 

17 

PE  TIJPN 

054  1  0 
054c'  OC 
054  3  OC 

1 6  0 

IF  'N*  lfc-1,161.164 

05440 

161 

DO  168  1=1.6 

05450 

ti  •  I  -  =  PZ<  1) 

05460 

D'I+10.  =  OZ < I > 

05470 

168 

CONTINUE 

05480 
0549 0C 

60  TO  163 

05500 

164 

DO  165  1=1.6 

055 1  0 

D<  I  '  =  PI  •  I  > 

055c' 0 

D  <  I  ♦ 1 0  >  =  01 • I» 

055  3  0 

165 

CONT 1  NtJE 

0554  0 

163 

CONTINUE 

05550 

Tl  =  85.0-X 

0556  0 

T8=T 1*T1 

05570 

p  =  D '6>*T8*D<5> 

05580 

DO  170  1=1.4 

05590 

J  =  5-1 

05600 

P  =  P*T8+D ' J' 

056 1  0 

170 

CONTINUE 

05680 

0  =  D < 1 6 ' ♦T8+D ■ 15' 

05630 

DO  171  1=1.4 

0564  0 

J  =  5-1 

05650 

0  =  a*Tc"*-D  •  J+l  O  ' 

05660 

171 

CONTINUE 

05670 
0568  0C 

0  =  C'*T  1 

05690 

T4  =DSQRT  <X«PI> 

05700 

T6  =  SIN  <y:> 

05710 
0578  0C 

T7  =  COS  <X) 

05730 
0574 0C 

IF  ■ N>  180.180.185 

05750 

180 

T5  =  r  .  p-i*'  <  *T6  +  cP+C 1 

05760 

GO  TO  99 

05770 

185 

T5  =  <<F*0>* T6  -  'F -O' 

05780 

99 

Y  =  T5 

05790 

RETURN 

05800 

END 

05810 

SUBROUTINE  PRRT 

0535 OC 
0583 OC 


►SUBROUTINE  PRRT  -  N-LRYER  ELRSTIC  SYSTEM 


05340 

COMMON 

RMCOY' RP  < 1 0> . 

22(10)  . 

E  <5>  . 

V  <5)  . 

HH  •  4 )  . 

05850 

i 

H  <4  •  . 

RZ  <396 > . 

R  <396. 5) 

.  B  <396.5) 

.  C <396.5) 

0586  0 

d 

D  <396*  5> • 

RJ  <396  > > 

RJ1 <396) 

.  R JO <396) 

.  TITLE <80 

05370 

% 

TEST 'll), 

BZ  <1  00>  . 

X  <5t 4«  4> 

,  SC  <4) . 

FM  <’£ » * 

05880 

& 

F'M  .  4.4.4  ' 

•  P. 

Z  • 

RR. 

NS. 

05390 

d 

N. 

L. 

ITN. 

RSZ. 

PSR. 

05900 

POM. 

PMU. 

SF. 

CSZ. 

C  ST, 

059 1  0 

:< 

C  SR. 

CTR. 

COM. 

CMU. 

PS  I* 

0598  0 

d 

NLINE. 

NQUTP . 

NTEST . 

I. 

ITN4. 

059 "■  0 

V- 

¥  • 

LC. 

JT. 

TZZ. 

PR, 

05940 

d 

PR. 

P. 

EP. 

TIP. 

TIM. 

0*^5  ft 

Tl. 

T2, 

T  ?« 

T4« 

T5« 

l  h 

T6« 

T8P. 

T8M. 

UH« 

B.U  . 

ii*.  ’■*  T  ft 

■  B  J  0  • 

ZF, 

.::i , 

'  Z8« 

'61 « 

n*v-4.x  ft 

'  68 « 

PH« 

»H8. 

VF  8* 

VF  P8 

VF  4 . 

VF  P4 , 

VFFA. 

RUT. 

rd: 

06  i.i  ii  i'i  [ih  Tm  ijj  0 . 84 1  1  36  51,  i\£  -  n .  33999 1  04 

fit.  i.i  1  fi  4  ZF  =  W 

fit  08  fi  NTE'T  =  2 

04030  IF  'Pi  St  St  9 

04040  9  CONTINUE 

04050  NTE'T  =  HP -  P  *  .0001 

04 04 0  IF  ’ NTEST>  6,6,5 

06  i.i 70  4  C  ONT  I NUE 

04030  NTEST  =  PAP  ♦  .  0001 

04090  ZF  =  P 

04100  5  CONTINUE 

04110  MTE I T  =  NTEST  +  1 

0414:0  IF  1 NTEST-1 0>  8,8,7 

04130  7  CONTINUE 

04140  NTEST  =  10 

04150  $  CONTINUE 

04 1 4 OC  ♦♦  COMPUTE  POINTS  FOP  LEGENDRE-6AUSS  INTEGRATION  ♦♦ 

04170  15  F  =  1 

04180  CALL  CHECK  <9> 

04190  ZF  =  2. 0*ZF 

04800  $Z2  *  0.0 

04210  00  28  1=1. I TN 

04820  $21  =  SZ2 

048  30  SZ2  =  B2  < I ♦ 1 > -ZF 

04240  2F  =  328  -  3Z1 

04850  PP  »  328  ♦  321 

04840  $G1=3F*61 

04870  3G2*3F*G2 

04890  AZ  f K i =PP-3G1 

04890  AZ<K>1:«  =PP-SG2 

04 3  0 0  AZ 1 F  *c > =PP>  3G8 

04 3 1 0  AZ '  K ♦ 3 > =PP+3G 1 

04  380  F  =  K  ♦  4 

04330  CALL  CHECK  (10) 

04340  88  CONTINUE 

04350  40  RETURN 

04360  END 

04370  SUBPOUT  I NE  CALC  IN  < I RT > 


04380C 
0639 OC 

♦♦♦♦♦♦$ 

UBROUT1NE  CALC  IN  - 

N-LAVER 

-♦ 

04400 

CDMMON 

'PMCOY^RRCIO), 

22  <  1  0>  * 

E  <5> ,  V • 5> , 

HH  f 4> , 

064 1  0 

& 

H  <4>  . 

AZ <396* , 

A <396, 5),  B  <396*  5> 

,  C  f 396, 5>  , 

0442  0 

* 

D<396«5> « 

A J  <396) , 

RJ1 <396> ,  RJ0<396> 

,  TITLE  <80> , 

044  30 

•J: 

TEST  <’  1 1>  , 

B2<100> , 

X  <5 » 4 1 4  > f  SC  *  4>  9 

FNf8,2>  , 

06440 

8. 

PM <4, 4, 4 > 

•  R  > 

Z  9  RRr 

NS, 

06450 

% 

N, 

L, 

I TN,  PS2, 

PSR, 

06460 

& 

RDM, 

RMU, 

SF,  CSZ, 

CST, 

06470 

%, 

CSR, 

CTR, 

CON,  CMU* 

PSI, 

06480 

8. 

NLINE, 

NOUTP* 

NTEST,  I, 

ITN4. 

Or.  4*  ft 

►  • 

LC  f 

JT,  TZ2, 

PR, 

06500 

V.' 

PA, 

p. 

EP,  TIP, 

TIM, 

065 1  0 

£v 

Tl, 

T8, 

T  3,  T4, 

T5, 

0658  0 

T6« 

T2P* 

T2M,  UA, 

B  J 1  • 

Ot  5  3  0 

'i 

B.IO, 

ZF, 

3Z1 ,  328, 

361, 

06540 

t 

3G8, 

PH, 

PH8,  VF8, 

VKP8, 

B-18 


0655  o 
0656  C> 
06570 
CiiJ5t.fi 
06-590 
06-6  0  Cif 
066 1  0 
066  c.' 0 

066  30 
06640 
06650 
06660 
06670 
06660 
066?  0 
06  7  0  0 
06  710 
06-720 
06720 
06  74  0 
06750 
06760 
06770 
067600 
06790 
06800 
068 1  0 
0688  0 
068  50 
0684  0 
0685  OC 
0686  0 
06870 
0688  0 
06890 
06?  0  0 
06?  1  0 
06980 
06930 
0694  0 
06950 
0696  0 
06970 
0698  0 
06990 
0  7  0  0  0 
07  0 1  0 
07020 
070  30 
0704  0 
0705  0 
07  060 
07070 
07080 
07090 
07100 
071  1  0 
07180C 
071  30 
07140 


¥  Vk.4  ,  VKP4.  VKK8.  POT. 

C OMMON  DPT  - C OUNTP . DEF , 4 > , NLP VS 
INTEGER  CDUNTP 
DIMENSION  l.l  ■  4  > 

DATA  1.1  0.  34785485- 2*0.  65214515.  0.  34785485 

2  VL=2. 0*V<L  > 

EL= ■ 1. 0+V  'L 1  ■ - E  ■  L  > 

VL 1=1. 0-VL 
C  22=0. 0 

c :t=o. 0 
csp=o.  0 

CTP=0. 0 
C OM=  0 . 0 

crin=o.  0 

NT  2  1  =  NTEST  ♦  1 

it:  =  1 

JT  =  0 

ARP  =  AR 

IF  - NOUTP  >  4.4.5 

4  ARP  =  ARP*PSI 

5  CONTINUE 

10  HO  40  1=1. ITN 

INITIALIZE  THE  SUB- INTEGRALS 

fcSZ=0.  0 

RST=0.  0 

RSR=0. 0 

PTP=  0. 0 

PON=0. 0 

RMU=0. 0 

COMPUTE  THE  SUB- INTEGRALS 
K  =  44(1-13 
DO  30  .1=1.4 

.11  =  1  ♦  J 
P=AZ ' J 1 • 

EP=E>:P  .p*Z> 

T 1  =B  (  J1 .  L  1  *£P 
T8=D ' J 1 . L 1 •  EP 
T  1F'=T  1  +T8 
T 1M=T 1 -T8 

T 1  =  ■  A  '  J  l«L3>B<Jl«L3  *Z>  *EP 
T2=  <C < J 1 . L> +D •  J 1 » L  • 42)  EP 
T2P=P* • T 1+T23 
T8M=P*<  T1-T23 
Wfi=H  J  I  J  1  '♦•»!(  J,» 

CALL  CHECK  <3> 

IF  (P>  80.80,15 
15  B  J 1 =RJ 1 •  Jl>  *P 
B  J0=P  JO  >  J  J  .3  *P 

PSZ=P$Z*WA*P*BJ0*<VL1*TIP-T2M> 

ROM=ROM+i,IA<>£L*B  JO*  <2.  0*VL  1  *T1M-T8P> 

PTP=PTR+WA*P*B J 1 ♦ i VL *T 1 M+T8P3 
PMU=PMU+WA*EL*BJl*  f T 1P*T£M> 

RSP=P$R*WA* . P*B JO* < < 1 . 0*VL>  *T1 P*T8M> -B J 1  * CT 1 P* T2M)  .-'P> 
PST=P<T*WA*i VL*P*BJ0*T1P*BJ14<T1P*T8«> -P> 

CALL  CHECK  <4* 

GO  TO  30 

SPECIAL  ROUTINE  FOP  P  =  ZERO 
80  PP=P*P 

RSZ=RSZ*WA*PP* i VL 1 *T 1 P-T2M- 
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M  7  1  *=■.  M 

60M=^0M*I.ih*EL*P»  <  3.  II ♦VI  1  ♦  T1 M  -  T8F  ■ 

07160 

PI  T  =PI  T ♦i.Ih#F'P#  •  1  VL  +  ii.5  •  •TIP* II.  5#T8M' 

1171  70 

pif-pit 

07 1  SO 

call  check  ■  5  ■ 

0  71?  0 

0  CONTINUE 

07c'  0  I’m; 

07  c' 1  0 

IF  =  •  AZ  ‘K.+4  •  -  AZ  •  K  ♦  1  >  >  1 . 7333736 

0  735'  0 

C $Z*CSZ+PSZ+iF 

07330 

C'TCST+PST^SF 

0734  0 

C  3  P = C  I P  *  P  S  P  ♦  S  F 

07350 

CTP=CTP+PTP*5F 

0736  0 

COM=COM+POM^SF 

0737  0 

CMU=CMU+RMU»SF 

0733  0 

PSZ  =  3. 0*PSZ»AP»SF 

07390 

CALL  CHECK  <6> 

07300 

TEITH  =  PBS  -PSZ>-10. 0*»<-4» 

07  3 1  0 

IF  •  ITS-NTS1  >  31.33.33 

07  330 

31 

CONTINUE 

07  3  30 

TEST -ITS'  =  TEITH 

07  34  0 

IT!  =  ITI  +  l 

07  350 

CO  TO  40 

07  36  0 

3«5 

CONTINUE 

07370 

te: T • NT i 1 i  =  TEITH 

07380 

DO  33  J  =  l.NTE'T 

07  390 

IF  • TEC  TH-TE IT  <  J)  >  35.  36.36 

07400 

-.er 

C ONTINUE 

0741  0 

TEITH  =  TE- T • J> 

07430 

CONTINUE 

074  30 

TE3T 1 J 1  =  TE'IT<J+1> 

07440 

CONTINUE 

07450 

IF  •'  TESTH  •  50.50.40 

07460 

40 

CONTINUE 

07470 

JT  =  1 

074:30 

CALL  HICHM 

07490 

50 

csz=c :z*app 

07500 

CALL  CHECK  -7' 

07510 

c :t=c IT  ♦APP 

07530 

C  TP«CTR4RPP 

07530 

C IP=C IP ♦AFP 

0754  0 

COM=COM*APP 

07550 

IF  "IOUNTP.  NE, 99>  PEF  < IPT>  =COM* 1 000. 

07560 

CMU=CMU*APP 

07570 

bit;  =  c iz+cst+c ip 

07580 

VI  TP  =  • C I Z- V • L  > ♦ *  C I  T +C  J  P> '  E 1 L ■ 

07590 

BIT  =  BITS  ♦  • 1 . 0-3. 0*V  <L ' >  E<L> 

07600 

IF  • TZZ '  73.78.71 

076 1  0 

n 

w  =  ~  7 

07630 

CONTINUE 

076  30 

phi  =  • c sp-v *  l • ♦ • c iz+c st> >  e  <l> 

07640 

-IT =8.  ('♦  <  1 . 0+V  ■  L  >  » ♦C  TP  E  *L  ’ 

07650 

PUT  =  OUT  -  V'L>  ♦  <CSZ  +  CIP'>.-E'L> 

07660 

IF  'COUNTP. NE. 99 »  CO  TO  99 

07670 

WRITE  '6.315'  P.Z.CSZ.CST.CMP.CTP.BSTS 

07680 

WRITE  C6»  318>  VSTP. PDT * PDS. SI T >  BIT 

07690 

WRITE  ■ 6.31 7>  COM 

07700 

31 

FORMAT  <1H0.F4. 1»F6. 1 » 1X.5HSTPE  . 1P5E11 

0771  0 

318 

FORMAT  • 1H  > 1 1X.5HSTRA  »1P5E11.3> 

07730 

317 

FORMAT  •  1 H  .  1 1 X .  5HDSPL  .  1  PE  1 1 . 3,' 

07730 

NLINE  =  NLINE  +  3 

07740 

IF  •  JT>  9*4 .  60 

H-20 


h  775 11  6- 11  W  P  I  T  E  '  6*  ,  3  1 6>  ■* 

or 760  316  FORMAT  •  t H+  .  1  37X.  4H SLOW < 

o r r r o  9*  return 

07780  END 

0  7  79  0  :  U  BF'OU  T  I  ME  HI  6HM 

073 OOC +  *  ♦♦♦♦♦♦SUBPOUT I HE  HI6HM  -  N-LAYEP  ELASTIC  SYSTEM  ♦♦♦♦♦♦ 

0781 OC  OCTAL  LIST 

0788 OC  LABEL 

07330CHIGHH  ♦♦♦♦♦♦SUBPOUT I HE  HU3HM  -  N-LAYEP  ELASTIC  SYSTEM  ♦♦♦♦♦♦ 

0784 OC 

07850  RETURN 

0  7860  EH  I' 

07870  SUBPOUT IHE  CHECK  <K IN> 

0788 OC  +  ♦  ♦♦♦♦♦♦SUBROUTINE  CHECK  -  N-LAYEP  ELASTIC  SYSTEM  ♦♦♦♦♦♦ 

0788 OC 


07800 

COMMON  ✓ PMCOY 

-  PR  '  1  0>  . 

ZZ • 1 0> . 

E  •  5>  . 

V  >  5  >  , 

HH '4 ■ , 

078 1  0 

H  '  4  >  » 

AZ  < 396  > , 

A <396, 5) 

,  B  1  3 86 ,  5  ^ 

•  L 1 5)  < 

07880 

D '386. 5} . 

A  J  *  386  ' , 

PJ1 ' 396 > 

,  RJO  <  396 ' 

,  TITLE <£0> 

07830 

y.. 

TEST  Cl 1>  » 

Bzaoo> , 

X  <5 , 4, 4> 

,  SC'4>, 

FM  ■  .£»£}  , 

07840 

V 

PM • 4*4, 4> 

» R, 

~r 
£-  J 

AR, 

NS, 

07850 

y. 

N» 

L, 

ITN, 

PSZ, 

RSR, 

07860 

% 

ROM. 

RMU, 

SF, 

csz. 

CST . 

07870 

t 

CSP, 

CTR, 

CDM, 

CMU, 

PS  I , 

07880 

NLINE, 

NOUTP, 

NTEST , 

I, 

ITH4. 

07880 

9:. 

K  , 

LC, 

JT, 

TZZ, 

PR, 

08  0  0  0 

y. 

PA. 

P, 

EP, 

TIP, 

TIM, 

08  0 1  0 

:t. 

Tl. 

T£, 

T  3, 

T4 , 

T5, 

08  08  0 

T6» 

T£P, 

T8M- 

WA , 

BJ1 , 

MX  M  -ill 

B  J  0 , 

ZE, 

321, 

SZ£. 

Sb  1  , 

Oo  04  0 

:< 

3G£» 

PH, 

PH£» 

VK£, 

VKPc , 

0£  nS  0 

VK4 . 

VKP4 , 

VKCK8, 

RDT , 

PBS 

flX  Or;.  M 
ijx  fjT  0 

MX  MX  0 

rx  08  0C 

RETURN 

END 

•JUFPDUTINE  cdfe 
u:f  fop  nll  i 

•K  IN' 

PROBLEMS, 

UP  TO  MAX 

DIMENSION 

OF  15  LAYERS 

081 OOC  PEPPOGPAMMEB  1  MAY  1880  BY  L  J  PAINTER  -  EXCELLENT  ACCURACY 


08 1  1  OC 

NOTE 

DOUBLE  ENTRIES  FOR 

C0E5  i: 

CO 15  PREVIOUSLY  USED 

08  IS  01 

08 1  3  OC 

•♦•♦♦♦SUBROUT INE  COEE  - 

N-LAYEP 

ELASTIC  SYSTEM  ♦♦♦♦♦♦ 

03 1 4  0 

COMMON 

PMCOY  RP  < I  O  ' , 

SZ(10)  , 

E  <5>  , 

V  <5>  , 

I 

X 

■U 

08  1 5  0 

V, 

H  <4  '  , 

AZ  <396} 

,  A  <396, 5> 

,  B <396,5}, 

C 1 396 • 

08 1 6  0 

:> 

D  1  386 ,  5 '  , 

A  J • 396> 

,  R.J1  <396) 

,  RJO' 396> , 

TITLE-, 

03 1  7  0 

TEST'- 1  1  '  . 

BZ  ■:  I  00' 

«  X  <5 , 4 , 4  } 

,  SC  <4  > , 

FM • 4} , 

08 1 8  0 

PM ■ 4  *  4 , 4 > 

,P, 

, 

AP, 

NS. 

08 1 8  0 

N- 

L, 

ITN, 

RSZ . 

RSR* 

ij  0 

POM. 

RMU, 

SF. 

CSZ, 

CST. 

0  85  1  n 

CSP. 

CTR, 

COM. 

CMU. 

PS  I . 

08550 

ML.  I  HE . 

NOUTP - 

NTEST. 

I , 

I TH4 « 

085  0 

y. 

K  . 

LC, 

JT. 

TZZ, 

PP. 

0854  0 

FA. 

P, 

EP, 

TIP. 

TIM. 

085*50 

Tl  - 

T£, 

T  3. 

T4, 

T5  > 

0828  0 

T6« 

T8P, 

T  £M« 

WA, 

BJ1 . 

0857  0 

B  JO, 

ZF. 

SZ1  > 

SZ£, 

SGI , 

U8  5  y  0 

y 

S  6£» 

PH, 

PH£, 

VK  £ , 

VK.P£, 

UXc'7  0 

08 0  0 

08 1 1  0 

08  ‘<5  0 

PERL04 

ENTF'V 

ENTRY 

vh  4 . 

0  •  £ ,  £  ■ 

COES '>  IN • 

C  0 1 5  *  K IN 

VK  P4 . 

VK  K  8  * 

RDT, 

RDS 

08330  LC  =  KIN 

08 34 OC  S -MX  SET  UP  MATRIX  X  =DI^MI^K I^K*M^D 


B-21 


rift  ar 

COMPUTE  THE  MmTPICE;  1- * p  ' 

ny  36  0 

1 

DO  10  K  =  1  >  M 

08  37  0 

T  l  =E  •  K  •  ♦  ‘  1  .  O+V  •  1  ♦  1 )  ■  <  E  -  K  *  1  *  *  ■  1 .  O+V  •  K  >  '■  1 

I'lft  -CM  II 

T1M=T1-1. 0 

lift  3ft  tl 

PH=P4H ' K  ■ 

054  0  0 

PH2=F'H*2.  0 

084 1  0 

VK 2=2. 0*V '  K  ■ 

084c'  0 

vi >  ps=a.  o*v  ■  y  *t  > 

054  3  0C 
0544  0 

VK K 8=8 . 04V  •  K > 4  V ■ K ♦ 1 > 

08450 

VKP4 =2. 0*VKP2 

0848  0 

VK 4=2 . 04VK2 

08470 

X ■  K  » 1. 1 > =VK4-3. 0-T1 

0848  0 

X  <K  j  2  <  1  >  =  0 .  0 

08490 

j  >  =T1M*  *:PH2-VK4+1 . 0> 

055  0  0 

<K  .  4-  1  ■  =-2.  04 T  1M*P 

055 1  OC 
08520 

T  3=PH24 ■ VK  2- 1 • 0» 

08580 

~4=VKK8* 1 . 0-3. 0*VKP2 

0854  0 

T5=PH2'4.  VK  P2- 1  .  O  ' 

08550 

T6=VK K  3+ 1 . 0- 3 . 04VK  2 

0858 or 
0357  0 

X  <K  -  1  •  2  '  =  ■  T  3+T4-T  1 4  « T5+T6  •  >  P 

08580 

X  '•  Y .  £ .  2  •  =  T  1  ♦  <  VK  P4-  3 .  0 >  - 1 .  0 

03590 

X  •  K  »  4*  2>  =T  1114  •  1  .  0-PH2-VKP4) 

086  0  OC 
086 1  0 

v , k , 3. 4  *  =  • T3-T4-T 1 4  <T5-T6  > >  P 

0862  OC 
086  30 

T  3 = PH24PH - VK K  8  + 1 . 0 

0884  0 

T 4*PH2'4  <  VK2'-VKP2> 

0885  OC 
08880 

X  <K  ,  1»4.>»  • :  T  3+T4+VKP2-T 1 4  <  T3  +  T4+VK  £>  )  'P 

08870 

X <K«  3. 2  >  =  <:-T3+T4-VKP£4Tl4<T3-T44VK£>>  P 

0868 OC 
0869  0 

X  <K  .1*3)  *T1M4  <  1 . 0-PH2-VK4) 

08700 

: :  <K ,  2 .  3) =2 . 04  T 1 M4P 

08710 

X<K. 3. 3>=VK4-3. O-Tl 

08720 

X<K*4.3>»0. 0 

08730C 

08740 

X  <K .  2. 4  •  =T  1114  CPH2-VK P4+ 1 .  O' 

08750 

X  CK  ,  4 . 4  >  =T  1 4  C  VK  P4-  3 .  0  >  - 1 .  0 

0876 OC 

K  =  K 

08770 

10 

CONTINUE 

0878 OC 

COMPUTE  THE  PRODUCT  MRTRlCEi  PM 

08790 

SC  <N) *4.  04. VIS'  -1.  O' 

08300 

IF  CM- 2 '  13- 11*11 

083 1  0 

1 1 

DO  12  K 1 =2- M 

03320 

M=NS-K 1 

08330 

C  C  <  M  *  =SC  1  M4 1 ' 44 . 4  .  V  •  M  >  -  1 .  ’ 

08340 

l£ 

C  OUT  I MIJE 

08850 

13 

COMTIMUE 

0886 OC 
08870 

C"  1  •  1  ■  =  1  . 

08880 

C  C2 . 2 » =  1 . 

08890 

0  c  1 , 2 '  =  0 . 

08900 

OQ  =P42.4H'N> 

08910 

IF  <00.  LT .  -88.  »  C'0=-88. 

03920 

IF  <00-88 .  ■  15-15-16 

08930 

15 

CONTINUE 

0894  0 

0  *'  1  •  2  1  1  —  00  * 
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AD-A099  593  ARMY  ENGINEER  WATERWAYS  EXPERIMENT  STATION  VICKSBURG— ETC  F/G  13/2 

NONDESTRUCTIVEE  TESTING  FOR  LIGHT  AIRCRAFT  PAVEMENTS.  PHASE  II.— ETC(U> 
NOV  80  A  J  BUSH  D0T-FA78WAI-898 

NL 


UNCLASSIFIED 


FAA-RD-80-9-2 


1.  i:  HOT  NEEDED  FOP  INITIALIZING  THE  PM  MATRIX 


Ij  6  ■*» 

16 

continue 

fi::'9?0C 

fl  LOOP  INITIALIZE;  PM-N-  .  • 

0898  fi 

DO  20  M= 1  « 4 

I'lS'J'j  fl 

LL=  '  M+  1  >  <- 2 

\y-*  0  0  0 

DO  20  J=3*4 

09  fl  1  fl 

PM  •  N*  M*  J  ■  =  X  •  N * M •  J  •  •  U  ■  L  L  •  2  • 

0’?*  0  8  0 

80 

CONTINUE 

09  1.1 1  0 

DU  26  K  1=2,  N 

09  i  >4  0 

Y =n:-f 1 

O  r*  05  0 

YY=Y+l 

0 1*06(1 

00  =P*2 .  ♦H  <K  > 

1*1=1 070 

I F • 00, LT. -88. >  00= -88. 

ITS  fi i'i 

I F  *■  00-88 .  .>  88  *  88  9  83 

O'?  09  0 

22 

CONTINUE 

09 1  o  o 

0  *"2*  1 '  =EXP  •" 00 j 

09 1  1  0 

0  <  1  *2>  =  1  -  -'P<2,  1> 

09 1 2  0 

60  TO  24 

09 .1  3  0 

23 

CONTINUE 

O'*  140 

P  •'  1  *  2>  =  0 . 

09 150 

Q*2* 1 >  =  1 . E20 

09160 

24 

CONTINUE 

09170 

DU  25  M= 1  *  4 

09 1 8  0 

LL=CM+1'  '2 

09190 

DU  25  J=3.4 

09200 

PM  'K  •  M*  .*>  =  <  X<K*M*1>  ♦  PM  •  KK *  I  *  J  • 

092 1  0 

v> 

♦X  CK  *  M*  2>  ♦  PMCKK,2,  J>  >  ♦  CKLL.  1> 

09220 

J? 

+  f  XOC*M*  3>  ♦  PM CKK,  3*  J> 

09230 

;< 

+X  CK 1  M *  4  ■  *  PM CKK* 4* J>  >  ♦  QCLL*2> 

0924  0 

25 

CONTINUE 

092*50 

£*8 

CONTINUE 

0926  OC 
092  70 

SOLVE  FOP  C  CNS>  AND  P • NS> 

0923  0 

T3=2. u*V • 1 > 

09  c' 90 

T4  = T  3  - 1 . 0 

09  3  0  0 
09  3  1  0 

FM<1’=  P*  ' PM  <1*1* 3 • ♦  PM C 1  *  3*  3) >  ♦  T3*  CPM < 1  * 2* 3> 

-PM  <1,4* 3> ) 

09320 

FM  c£)  =  P*  • PM*  1  .  1  *  3>  -PMC1  *3*  3>  )  ♦  T4*  *  PM  0*2,3)  4PM  <1*4*  3) 

09830 

FM  •  3>  =  P*  <PM  0,1*4)  ♦PM  Cl  ,  3*  4>  )  ♦  T3*  <PM  ■  1 » 2*  4> 

-PM  Cl. 4* 4)) 

09  34  0 

FM  <  4>  =  P*  < PM  O  *  1*4)  —PM  O  *  3*  4>  )  ♦  T44  CPM  (1  *  2*  4)  +PH  <1  *  4*  4 >  > 

0935  0 

DFAC=SC  <  1  >  ✓  <  <FM  ■:  1  >  *FM  <4>  'FM  <3>  »FM  <2>  )  ♦P^P) 

09  36  0 

A  CLC  *  NS)  =  0.0 

09370 

B  CLC  •  NS>  =  0.0 

09  38  0 

CvLC*NS)  =  -PM '3' ♦DP AC 

09390 

D  <LC  » NS>  =  FM  Cl > ♦DFAC 

094 00C 

BACK  SDL VE  FOP  THE  OTHEP  A*B.C.D 

094 1  0 

DO  91  K  1  =  1 ,  N 

0942  0 

A  CLC  j  K  1 )  =  <PM  CK  1*1*  3)  *C  <LC  *  NS>  ♦PM  CK  1*1*  4>  ♦!'  CLC  * 

NS)  )  '■  SC  CK1  > 

09430 

B  CL  C  *  K 1 ) * CPM  CK  1  >  2 ,  3) *C  CLC  *  NS) +PM CK 1  *  2 « 4> ♦D  <  LC  * 

NS)  )  ''SC  CK1) 

1-19440 

C CLC »K1>» CPM  CK 1*3* 3> *C CLC  *  NS)  +PM  CK 1  *  3  *  4) ♦D<LC, 

NS) ) ✓SC CK  1 ) 

09450 

91 

D  CLC  *  K 1 >  =  CPM  CK  1 , 4 .  3  '  ♦C  CLC  *  NS)  ♦PM  CK  1  *  4 , 4>  *D  CLC  * 

NS) > ✓SC CK 1 ) 

09460 

1  00 

CONTINUE 

094  70 

PETURN 

09480 

END 

09490 

SUBROUTINE  SIMEQ  ■ A* B* N* KERP* IDIM) 

095  0  OC 

SIMEQ 

095 1  OC 

THIS  SUBROUTINE  SOLVES  A  SYSTEM  OF  LINEAR  EQUATIONS 

0952 OC 

AX=B  BY  THE  METHOD  OF  PIVOTAL  CONDENSATION. IT 

IS  USED 

095  3  OC 

FOP  DENSE  MATRICES  OF  COEFFICIENTS. 

0954  OC 

CALLING  ARGUMENTS: 

0'S  5  5  fir 

A:  THE  NAME  OF  AN  N  BY  N  MATRIX  OF  COEFFICIENTS  OF  THE  EQUATIONS 

0958 OC 

WHICH  IS  DESTROYED  DUPING  COMPUTATION 

. 

B-23 


0857 OC:  f:  THE  NAME  OF  AN  ARRAY  CONTAINING  THE  N  CONSTANT; . 

Mi; :  H:  THE  ORDER  OF  THE  SYSTEM 

0858OC:  CEPR:  INDICATOR  RETURNED  BY  THE  SUBROUTINE  WHICH  II 

iWhimiI:  ONE  IF  :  VS  TEH  IS  SINGULAR  AND  ZERO  OTHERWISE. 

o?6  1  *;*•: :  feaii  a  column  by  column. 

O'?*:  8  oc  : 

086  5  <>  dimen:  ION  h-  I  DIM,  IDIH‘  .F  -  IPIM>  .L  <50*  ,M«50> 

o?f-4M  eps=ie-6 

m-=i65o  ep:op=fp:*ep: 

iY*66  M  t  EPR=0 

086 7 OC:  CLEAR  OUT  FERMUTAT I  ON  VECTORS 

0*680  DO  3  J  =  1  •  N 

0*6  8  O  M  •  I  *  =  i  t 

Cr?  700  5  L  •  I  •  =  0 

0571 OC !  LOOP  FOF  N  PIVOT  POINTS 

0*780  DO  14  F  P=  1 « M 

iV?7  5  0  F  =  0. 

fr?740  F  ;C'F=0. 

O'?  7?  0  DO  7  I  =  1  •  N 

0*760  DO  7  ,i=l  •  N 

Cr?77  0  I F  •  M  1  I  ■  >7*4.7 

0878  0  4  I F  i  L  •  J  •  ■  7 •  5 «  7 

08780  S'  T=A'l«.i' 

n?:500  T  ;op=T*T 

O'?;? i o  if  •t:op-p:op>?»7»6 


i.i -f* Sc  0  £ 

P  =  T 

1 1  *•:  1 1 

p : cr= t : op 

0'?:54  0  1P=I 

ivixs.  0  ►  C  =  J 

0886 0  7  CONTINUE 

0*870  IF  ■  PSC'P-EPSC'R'  17.  1 7.  8 

08880  8  h  1  ►  P  •  =('  C 

0*880  L>KC'=1 

088 OOCs  DIVIDE  Y EV  POM  BY  PIVDT 

08810  DO  30  J= 1 • N 

0888 0  i p  <L ■ J ■ > 3  0 >  8  >  3 0 

088 3 0  *  A • K  P .  J • =A • ► P « J »  -  P 

09940  90  CONTINUE 

088?  0  BKP=B‘t,P»'P 

0888  0  B 1 K  P  >  =  BKP 

0887 OC  S  SUBTRACT  MULTIPLE  OF  KEY  POW  FROM  OTHER  ROWS 

08880  DO  14  1*1  *N 

08880  IF  I-KP*  31*  14.31 

10000  31  P=A • I « KC > 

10010  DO  33  J  = 1 <  N 

10080  IF  <L* J*»  33. 38.33 

l 0  0 3  0  38  A ■ 1 1 J • =A ' I . J > -R*A CKP . J > 

1004U  33  CONTINUE 

I  0 0? 0  B  •  I  '  *B  •'  I)  -P*BKR 

10060  14  CONTINUE 

10070C*  REORDER  RESULTS 

10080  DO  3?  I  =  1 .  N 

10X'  ?0  I  P=M  I  • 

10100  3?  A > IP. 1 >  =B '  P 

10110  DO  36  1=1. N 

10180  36  B  >  I  >  =  A  •  1 . 1  > 

10130  RETURN 

10140CS  ERROR  ACTION 

10150  17  KERR* 1 

10160  RETURN 

10170  END 


APPENDIX  0:  GUIDE  TO  USE  OF  COMPUTER  PROGRAM  PAVFVAL 

The  computer  program  PAVEVAL  calcul ate::  I. lie  allowable  load- 

carrying  rapacity  and  the  required  overlay  thickness  for  rigid  and 

flexible  pavements.  A  program  listing  is  contained  in  the  report  by 
3 

Weiss.  Input  guides,  typical  inputs,  and  typical  outputs  are  fur¬ 
nished  in  this  appendix  for  evaluation  of  load-carrying  capability  of 
both  flexible  and  rigid  pavements. 

INPUT  GUIDE  FOR  FLEXIBLE  PAVEMENTS 

Line  1 

Title  80  characters 


Line  2 


NSYS 


NSYS  =  Number  of  problems  to  run 


Line  3 


EKEY 


EKEY2 


EKEY  =  Limiting  strain  and  stress  subroutine  code:  set  =  2, 
for  calls  to  subroutine  flex 

EKEY2  =  Pavement  problem  code:  set  =  0,  for  allowable  load 


Line  U 


ES 

EA 

YRN 

ALOAD 

ALIN 

CAREA 

DSM 

SWL 

PCRATIO 

ES  =  Subgrade  modulus,  psi 

EA  =  Asphalt  modulus  of  existing  layer 

YRN  =  Yearly  load  repetition  number 

ALOAD  =  Initial  load,  lb 

ALIN  =  Load  increment,  lb 

2  2 

CAREA  =  Contact  area  (irr  ),  in. 


C-l 


DSM  =  Dynamic  stiffness  modulus,  for  reference  only  (any  number) 
SWL  =  Set  *  0 

PCRATIO  =  Pass-to-coverage  ratio 


Line  5 


NLAYS 


ISMO 


IRED 


NLAYS  =  Number  of  layers  in  payement  system 
ISMO  =  Request  for  rough  computational  procedure:  set  =  0 
IRED  =  Input  format  for  Line  6:  set  =  0 
Line  6-1  through  6-(NLAYS-l)  (one  for  each  layer,  except  last  layer) 


E(i) 

NU(  i ) 

THICK( i ) 

AK(i) 

E(i)  =  Modulus  of  layer  i 

NU(i)  =  Poisson's  ratio  of  layer  i 

THICK(i)  =  Thickness  of  layer  i 

AK( i )  =  Interface  compliance:  set  all  AK(i)'s  =  0 


Line  7 


E( NLAYS) 


NU( NLAYS 


E( NLAYS )  =  Modulus  of  last  layer 
NU(NLAYS)  =  Poisson’s  ratio  of  last  layer 


Line  8 


NLOAD 


NLOAD  =  Number  of  loaded  areas 
Line  9  (one  for  each  load) 


LOSTRS(i) 

|  RADIUS (i ) 

X(i) 

Y(i) 

HOSTR(i) 

PSI(i) 

LOSTRS(i)  =  Vertical  load  for  load  area  i 
RADIUS(i)  =  Radius  of  loaded  area  i 


C-2 


X(i)  =  Abscissa  of  center  of  loaded  area:  set  =  0 

Y(i)  =  Ordinate  of  center  of  loaded  area:  set  =  0 

HOSTR(l)  *  Horizontal  load  for  load  area  i:  normally  =  0 

PSl(i)  =  Angle  of  HOSTR(i)  with  respect  to  positive  X-axis  in 
in  degrees:  normally  =  0 

Line  10 
NPOS 


NPOS  =  Number  of  depths  that  will  be  used  for  iteration  purposes: 
set  =  2 

Line  11-1  through  11-NP0S 


LAYER ( I ) 

AX(  i ) 

AY(i) 

DEPTH ( i ) 

ETA(i) 

LAYER(i)  =  Layer  number  for  position  i:  set;  LAYER(l)  =  1,  Layer(l) 
=  NLAYS(last  layer) 

AX(i)  =  Abscissa  of  position:  set  AX(l)  =  AX(2)  =  0 

AY( i )  *  Ordinate  of  position:  set  AY(l)  =  AY(2)  =  0 

DEPTH(i)  *  Depth  from  pavement  surface  to  position:  set;  DEPTH(l) 

-  THICK(l),  DEPTH(2)  =  distance  from  pavement  surface 
to  top  of  subgrade 


ETA(i)  =  Angle  from  which  position  is  observed:  set  =  0 


TYPICAL  INPUT  FOR  FLEXIBLE  PAVEMENT  EVALUATION 

01  0»8H 

02038  I DENT s  POt  F  300?  EU3H 
0  •:  01!  OPTION:  FCIPTERN 
04 O  F  8  M ' S3  * 06  043  0/ 1 7 0 0 
0T.01iU.iE:  .STL  IT 
0>:038  FORTY:  XPEF 
0  V  038  SELECTIi:  F'PVEVRL 
03  03 :  MOTE :  F'ftVEVFiL 
090.3:  EXECUTE 

1 003*  L  MI  TS  8  30? 3 OK* . 15000 


110PTPF  1ft  SO  KIP  SINGLE 

LINE 

1 

ISO 

1 

LINE 

2 

130 

S  0 

LINE 

3 

140 

8320  7 70000 

25000  20000  1000  127 

1000  0  7.94  LINE 

4 

150 

3  0  0 

LINE 

5 

160 

770000  0.35 

7.5  0 

LINE 

6-1 

170 

40304  0.35 

SO. 5  0 

LINE 

6-2 

130 

8320  0.40 

LINE 

7 

190 

1 

LINE 

8 

soo 

9500  6.36  0 

0  0  0 

LINE 

9 

S10 

S 

LINE 

10 

sso 

I  0  0  7.5  0 

LINE 

11-1 

S 3  0 

3  0  0  S3 .  0 

LINE 

11-2 

24  03 :  END  JOE 
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INPUT  GUIDE  FOR  RIGID  PAVEMENTS 


Line  1 


Title  Bo  Oharar Lerc 


Line  2 


NY3YS 


NSYS  =  number  of  problems  to  run 
Line  3 


EKEY 


EKEY2 


EKEY  =  Limiting  strain  and  stress  subroutine  code:  set  -  3,  for 
calls  to  subroutine  RPAL 

EKEY2  =  Pavement  problem  code:  set  =  0,  for  allowable  load 
Line  L 


DSM 

FAC 

YRN 

R 

ALOAD 

ALIN 

CAREA 

SWL 


DSM  =  Dynamic  stiffnes  modulus,  for  reference 

FAC  =  Pass-to-coverage  ratio 

YRN  =  Yearly  load  repetition  number 

R  =  Flexural  strength,  psi 

ALOAD  =  Initial  load,  lb 

ALIN  =  Load  increment,  lb 

2  2 

CAREA  -  Contact  area  ( nr  ) ,  in. 


SWL  = 
Line  5 


set  =  0 


NLAYS 


ISMO 


IRED 


NLAYS  =  Number  of  layers  in  pavement  system 
ISMO  =  Request  for  smooth  computational  procedure: 
IRED  =  Input  format  for  Line  6:  set  =  ] 

C-1Y 


set  =  1 


Line  6-1  through  6-(NLAYS-l)  (one  for  each  layer,  except  last  layer) 


E(i) 

NU  ( i ) 

THXCK( i ) 

ALK(i) 

E(i)  =  Modulus  of  layer  i 

NU(i)  =  Poisson's  ratio  of  layer  i 

THICK(i)  =  Thickness  of  layer  i 

ALK(i)  =  Reduced  interface  compliance:  set;  ALK(l)  =  1000,  all 
other  ALK(i)'s  =  0 


Line  7 


E(NLAYS) 


NU(NLAYS) 


E(NLAYS)  =  Modulus  of  last  layer 
NU(NLAYS)  =  Poisson's  ratio  of  last  layer 


Line  8 


NLOAD 


NLQAD  =  Number  of  loaded  areas 
Line  9  (one  for  each  load) 


LOSTRS(i) 

RADIUS ( i ) 

X(i) 

Y(i) 

HOSTR(i) 

PSI(i) 

LOSTRS(i)  =  Vertical  load  for  load  area  i 

RADIUS(i)  =  Radius  of  loaded  area  i 

X(i)  =  Abscissa  of  center  of  loaded  area:  set  =  0 

HOSTR(i)  =  Horizontal  load  for  load  area  i:  normally  =  0 

PSl(i)  =  Angle  of  Hostr(i)  with  respect  to  positive  x-axis  in 
degrees:  normally  =  0 


Line  10 


NPOS 


NPOS  =  Number  of  depths  that  will  be  used  for  iteration  purposes: 
set  =  1 


C-18 


Line  11-1  through  11-NPOS 


LAYER ( i ) 

AX(  i  ) 

AY  (  i  ) 

DEPTH ( i ) 

ETA ( i  ) 

LAYER( l )  =  I,.-iy«>r  number  i'or  position  i;  net.  =  1 
AX(l)  =  Abscissa  of  position:  set  =  0 
AY(l)  =  Ordinate  of  position:  set  =  0 

DEPTH(l)  =  Depth  from  pavement  surface  to  position;  set  =  THICK(l) 
ETA(l)  =  Angle  from  which  position  is  observed:  set  =  0 
TYPICAL  INPUT  FOR  RIGID  PAVEMENT  EVALUATION 

010«:*N 

03  01 :  I  DENT :  P  0SF3 0  0  •  BU i  H 
0  3  01 : OPT I ON : FDPTPhn 
04  0i s  M  i  i?3 :  Ob  048 0  -  1  7 0  0 


O'?  01 : 

USE: . GTLTT 

0  b  01  ■ 

FOPTY::  f  r.F 

07  01: 

SELECT*: PRVEVRL 

0301: 

NOTE: PR''  *'  v‘!4L 

O'?  01 : 

EXECUTE 

1  0  01 : 

LIMITS:  3  0  *  3  OF  *  , 

15000 

1 1  0PC 

C  POhH  SECTION 

o 

LINE 

1 

1  3  0 

1 

LINE 

2 

1  3  0 

3  0 

LINE 

3 

140 

1000  7.34  35000 

799  1 5  0  0  0  1  0  0  fi  1  £ 7  0 

LINE 

4 

150 

3  1  1 

LINE 

5 

1 6  0 

4000000  0. 15  6 

1  0  0  0 

LINE 

6-1 

170 

7378?  0.30  10 

0 

LINE 

6-2 

130 

13378  0.40 

LINE 

7 

1  ■?  0 

1 

LINE 

8 

3  0  0 

850  0  8.  36  0  0  0 

0 

LINE 

9 

3 1  0 

1 

LINE 

10 

33  0 

1  0  0  6  0  - 

LINE 

11-1 

3  3  01: 

END  JOE 

F»Ei=tr'Y 


0-19 
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